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C^ : ABSTRACT 

We present optical and near-infrared photometry and spectroscopy of the type la SN 1998bu in the 
Leo I Group galaxy M96 (NGC 3368). The data set consists of 356 photometric measurements and 29 
spectra of SN 1998bu between UT 1998 May 11 and July 15. The well-sampled light curve indicates the 
supernova reached maximum light in B on UT 1998 May 19.3 (JD 2450952. 8±0. 8) with B = 12.22 ± 0.03 
and V = 11.88 ± 0.02. Application of a revised version of the Multicolor Light Curve Shape (MLCS) 
method yields an extinction toward the supernova oi Ay = 0.94 ±0.15 mag, and indicates the supernova 
was of average luminosity compared to other normal type la supernovae. Using the HST Cepheid 
distance modulus to M96 (Tanvir et al. 1995) and the MLCS fit parameters for the supernova, we derive 
an extinction-corrected absolute magnitude for SN 1998bu at maximum. My = —19.42 ± 0.22. Our 
independent results for this supernova are consistent with those of Suntzeff et al. (1999). Combining 
SN 1998bu with three other well-observed local calibrators and 42 supernovae in the Hubble flow yields 
a Hubble constant, Hq — 64tg km s~^ Mpc~^, where the error estimate incorporates possible sources of 
systematic uncertainty including the calibration of the Cepheid period-luminosity relation, the metallicity 
dependence of the Cepheid distance scale, and the distance to the LMC. 

Subject headings: distance scale — supernovae: general — supernovae: individual (SN 1998bu) 
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1. INTRODUCTION 

Type la supernovae (SNc la) have recently assumed an 
elite status at the top rung of the cosmic distance lad- 
der. Comprehensive studies of SNe la have established 
their general spectroscopic and photometric homogeneity, 
along with quantifiable diversity. SNe la make reason- 
ably good standard candles (Kowal 1968; Sandage & Tam- 
mann 1993; Branch & Miller 1993) and excellent calibrated 
candles (Phillips 1993; Hamuy et al. 1995, 1996b; Riess, 
Press, & Kirshner 1995a, 1996a; Tripp 1998; Phillips et al. 
1999) which, combined with their high intrinsic luminosity, 
makes them superb indicators of very large distances and a 
powerful tool for cosmology (Branch 1998, and references 
therein). 

The downside is that SNe la are rare and fleeting, mak- 
ing them more challenging to study than other astronom- 
ical distance indicators. Nevertheless, recent applications 
of these "standard bombs" have been numerous. Within a 
few hundred Mpc, they have been used to measure the ex- 
pansion of the Universe: with recession velocities of their 
host galaxies and distances to individual SNe la, the Hub- 
ble law has been tested to redshifts z ~ 0.1 at high preci- 
sion (Hamuy et al. 1996b; Riess, Press, & Kirshner 1996a; 
Tammann 1998). SNe la have also been used to measure 
peculiar motions of galaxies and large scale flows (Tam- 
mann & Leibundgut 1990; Miller & Branch 1992; Jerjen 
& Tammann 1993; Riess, Press, & Kirshner 1995b; Riess 
et al. 1997a; Zehavi et al. 1998; Tammann 1998) as well 
as to provide constraints on the properties of extragalactic 
dust (Riess, Press, & Kirshner 1996b). At larger distances, 
SNe la serve as standard clocks: cosmological time dila- 
tion has been demonstrated by comparing light curves of 
distant supernovae with light curves from those nearby 
(Leibundgut et al. 1996; Goldhaber et al. 1997), as well 
as from the relative rates of spectral evolution (Riess et al. 
1997b). Most recently, SNe la have been used to measure 
luminosity distances at high redshift (0.3 < z < 1.0), from 
which the geometry and expansion history of the Universe 
can be determined (N0rgaard-Nielsen et al. 1989; Perlmut- 
ter et al. 1995; Schmidt et al. 1998). These measurements 
have likely tolled the death knell of standard (SIm — 1) 
cold dark matter (CDM) cosmology, strongly disfavoring 
the possibility of enough gravitating matter to flatten the 
Universe (Perlmutter et al. 1998; Garnavich et al. 1998a). 
Even more surprisingly, preliminary indications show that 
the expansion is accelerating at the current epoch (Riess 
et al. 1998a; Perlmutter et al. 1999), inconsistent with 
the idea that ordinary matter is the dominant component 
of the Universe's energy density. Alternatives such as a 
cosmological constant, or other forms of energy density 
similar to it, seem to be favored (Garnavich et al. 1998b; 
Perlmutter et al. 1999). Combined with constraints from 
the cosmic microwave background anisotropy power spec- 
trum, the supernovae have also given preliminary indica- 
tions for a Universe with zero global curvature (Garnavich 
et al. 1998b; White 1998; Lineweaver 1998; Tegmark 1999; 
Efstathiou et al. 1999; Perlmutter, Turner, & White 1999; 
Roos & Harun-or-Rashid 1999). 

These studies rely on SNe la only as bright, precise, rel- 
ative distance indicators. However, the determination of 
the expansion rate of the Universe, the Hubble constant, 
and the age of the Universe require an accurate absolute 



distance scale. Measuring absolute distances to SNe la re- 
quires calibrating them through objects on lower rungs of 
the distance ladder. The best distance indicator for this 
remains what it has been since the days of Hubble him- 
self: 5 Cephei variable stars. Cepheids are bright enough 
to be studied in nearby galaxies, including those which 
have hosted SNe la. In addition they can be found in sig- 
nificant numbers, which allows for a precise mean period- 
luminosity (PL) relation to be constructed and compared 
to nearby samples. In general Cepheids in external galax- 
ies are compared to Cepheids in the Large Magellanic 
Cloud (LMC), whose distance is calibrated through a vari- 
ety o f means (not always with the same result, cf. Section 
5.2.3| ). The importance of Cepheids to calibrate secondary 
distance indicators (including SNe la) prompted the cre- 
ation of the Hubble Space Telescope (HST) Key Project 
on the Extragalactic Distance Scale, which has measured 
Cepheid distances to a number of nearby galaxies (Freed- 
man et al. 1994, 1998 and references therein). Another 
HST effort has been underway to measure Cepheid dis- 
tances to galaxies that have hosted SNe la, to calibrate 
them specifically (Sandage et al. 1992; Saha et al. 1997, 
and references therein). At present the use of SNe la to 
measure the Hubble constant is limited by the paucity of 
local calibrators, not by the number or precision of dis- 
tances to objects in the Hubble flow. 

Nature has been kind by providing us with SN 1998bu 
in the Leo I Group galaxy M96 (NGC 3368), for which 
a Cepheid distance had already been obtained with HST 
by a third group, Tanvir et al. (1995). In this paper we 
calibrate the absolute magnitude of SN 1998bu through 
extensive optical and near-infrared photometry and spec- 
troscopy. In §2 we describe our observations and reduction 
procedure. In §3 we derive parameters of SN 1998bu, in- 
cluding the extinction along the line of sight as well as 
a quantitative estimate of the intrinsic luminosity of SN 
1998bu compared to other SNe la. In §4 we combine the 
properties of SN 1998bu with other Cepheid-calibrated 
SNe la and SNe la in the Hubble flow to determine the 
Hubble constant (i?o) and its statistical uncertainty. We 
discuss our results, sources of systematic error, and im- 
plications for the age of the Universe, tg, in §5. Finally, 
we conclude and summarize in §6. Independent observa- 
tions and analysis of SN 1998bu have been performed by 
Suntzeff et al. (1999), and we compare our analysis and 
results with theirs throughout the paper. Infrared and op- 
tical spectra and uncalibrated light curves of SN 1998bu 
have also been presented by Meikle & Hernandez (1999). 

2. OBSERVATIONS AND ANALYSIS 

2.1. Discovery 

SN 1998bu was discovered by the amateur astronomer 
M. ViUi on UT 1998 May 9.9 on CCD images of M96 (ViUi 
1998). The supernova was located at a = 10'^46'"46!03, 
S = +ll°50'07'.'l (equinox 2000.0), about one arcminute 
north of the host galaxy nucleus (Nakano & Aoki 1998). 
At the Center for Astrophysics we monitor new, bright 
supernovae spectroscopically with the Center for Astro- 
physics F. L. Whipple Observatory (FLWO) 1.5-m Till- 
inghast reflector and FAST spectrograph (Fabricant et al. 
1998), and photometrically with the FLWO 1.2-m reflector 
in optical and near-infrared passbands. Our photometric 



observations of SN 1998bu with the 1.2-m began on May 
11.1, with the new STELIRCAM near-infrared InSb array 
detector. The discovery of SN 1998bu ocurred during lunar 
bright time, and so the l.S-rn telescope was not equipped 
with the FAST instrument until May 15. 

High resolution spectra taken by Munari et al. (1998) 
with the Asiago Observatory 1.8-m telescope showed inter- 
stellar Na I D absorption from our Galaxy, as well as from 
M96 at a heliocentric radial velocity of 744.8 ± 0.3 km s^^. 
Low-disperson spectra of SN 1998bu were taken by Meikle 
et al. (1998) on May 12.9 and Ayani, Nakatani & Yamaoka 
(1998) on May 14.5, which revealed the supernova to be 
of type la about a week before maximum light. Our first 
spectroscopic observations of SN 1998bu with the FLWO 
1.5-m were taken on May 16.1. 

A type la supernova in a galaxy whose Cepheid dis- 
tance had already been measured by HST (Tanvir et al. 
1995) provided a unique opportunity; it was the opposite 
of the usual case, in which an HST Cepheid distance to 
a galaxy is measured specifically because the galaxy was 
an SN la host. Well-measured light curves are the key to 
SN la distances, so we undertook extensive photometric 
observations in the UBVRIJHK passbands. 

2.2. Optical Photometry 

The plurality of our optical photometric observations 
(29 nights) was obtained with the FLWO 1.2-m telescope 
-I- "4Shooter" CCD mosaic camera (Szentgyorgyi et al. 
1999). The 4Shooter consists of a 2x2 array of thinned, 
back-side illuminated, anti-reflective coated Loral 2048^ 
CCD detectors, situated at the f/8 Cassegrain focus. The 
pixel size is 15 /im, yielding a scale of 0'.'335 per pixel at 
the focal plane and a field of view of approximately llf4 
on a side for each chip, with total sky coverage of 0.15 
deg^. Our observations were taken in a 2x2 binned mode, 
so that the resulting images were sampled at 0'.'67 per 
pixel, well matched to the typical seeing (l'.'5-2" FWHM) 
achieved at this telescope. All observations of SN 1998bu 
with the mosaic were made on the same CCD (chip 1), 
which has the best combination of cosmetic characteristics 
and quantum efficiency. Two nights of observations were 
obtained on this telescope with the "AndyCam" instru- 
ment, a CCD camera with a single CCD, very similar to 
those that make up the 4Shooter. Both instruments have 
good ultraviolet and near-infrared response, which enabled 
us to make observations in the Johnson UBV and Kron- 
Cousins RI bandpasses. Our optical filters are constructed 
from Schott glass components, as recommended by Bessell 
(1990) for coated CCDs. The FLWO BVRI filter prescrip- 
tions are described by Riess et al. (1999); in general they 
match well the prescriptions of Bessell (1990), though the 
FLWO / filter extends to somewhat longer wavelengths. 
The FLWO U filter transmission is also a good match to 
the Bessell (1990) UX specification. 

Our FLWO 1.2-m observations of SN 1998bu are part of 
an ongoing supernova monitoring program at CfA. Super- 
nova observations are not well suited to the scheduled time 
allocation procedure that is typical at most telescopes. To 
follow supernovae, the time allocation committee autho- 
rizes us to enlist the generous aid of the scheduled ob- 
servers, asking them to devote a small fraction of observ- 
ing time (usually limited to ~ 20 minutes per night) to 
the SN program. We complement the monitoring observa- 



tions with scheduled nights (usually one night per month) 
to measure fainter objects and perform photometric cali- 
brations. We have been quite successful observing in this 
mode; a set of 22 SNe la light curves garnered as a result 
of this program has been presented by Riess et al. (1999). 

The FLWO 1.2-m is equipped with an infrared instru- 
ment during bright time, which provides useful IR super- 
nova data (obtained in a similar observing mode), at the 
price of bright-time gaps in our optical light curves. For 
this object we made a special effort to minimize these gaps 
by inviting observers at other institutions to participate. 
Six epochs of optical photometric observations in UBVRI 
were taken at the Michigan-Dartmouth-MIT (MDM) Ob- 
servatory 2.4-m Hiltner telescope and direct imager, with 
the thinned, back-illuminated, 1024^ "Charlotte" CCD de- 
tector, located at the f/7.5 Cassegrain focus and providing 
a 4f7 X A'.l field-of-view at 0'.'28 per pixel. Further ob- 
servations on two nights were obtained at the Kitt Peak 
National Observatory (KPNO) 0.9-m telescope, with the 
T2KA (20482) CCD detector at the f/7.5 Cassegrain focus, 
yielding a 23' x 23' field-of-view at 0'.'68 per pixel. Target 
of opportunity observations were also carried out on five 
nights during NOAO time at the WIYN Observatory 3.5- 
m telescope with the S2KB (2048^) CCD Imager at f/6.3 
with 0'.'20 per pixel and a 6f8 x 6.'8 field-of-view. 

Other sites also observed SN 1998bu and we report those 
data as well. We include CCD data from the Whitin 
Observatory 0.6-m telescope at Wellesley College (eight 
nights using a 1024^ CCD at f/13.5, with a scale of 
0'.'91 per pixel), the Gettysburg College Observatory 0.4- 
m (twelve nights using a 1024^ CCD at f/11, 0'.'84 per 
pixel) and the 0.76-m Katzman Automatic Imaging Tele- 
scope (KAIT) at Lick Observatory run by the University 
of California, Berkeley (fifteen nights with a 512^ CCD 
at f/8. 2, 0'.'63 per pixel). The detectors used in these ob- 
servations were not as blue-sensitive as the others, so only 
BVRI images were taken. Our observations of SN 1998bu 
continued until July 2 when it was too close to the setting 
sun to provide good photometric data. In total our optical 
photometric data set consists of 327 measurements of SN 
1998bu. 

All CCD observations were reduced (uniformly, begin- 
ning with the raw data) in the standard fashion, with bias 
subtraction, dark current subtraction (not necessary in 
most cases) and flat-field correction using the IRAF CCD- 
PROC package. Most of the observations were taken in 
non-photometric conditions, so we have performed differ- 
ential photometry with a sequence of six comparison stars 
in the supernova field, shown in Fig. ^. Comparison stars 
1 and 2 were calibrated on four photometric nights (two 
from the FLWO 1.2-m, one each from the MDM 2.4-m and 
Kitt Peak 0.9-m). Stars 3, 4 and 5 were not in the MDM 
field of view, so these were calibrated from three nights. 
To calibrate the comparison stars into a local standard 
star sequence, Landolt (1992) standard fields providing 
stars in a wide range of color were observed in UBVRI 
over a wide range in airmass. The supernova field was 
also observed in these filters at an airmass within the air- 
mass range of the standard star observations. These data 
were reduced and stellar instrumental magnitudes were de- 
termined from aperture photometry using the APPHOT 
package in IRAF. We then derived zero points and trans- 
formation coefficients linear in airmass and color from the 



standard stars using the prescription of Harris, Fitzgerald, 
& Reed (1981). This transformation was then appUed to 
the comparison star instrumental magnitudes to determine 
their standard magnitudes. We treated observations from 
each photometric night independently, and averaged the 
final standard magnitudes. Table |l| displays these mean 
standard magnitudes, along with the error in the mean 
determined from the scatter of repeat observations. We 
also list the comparison star identification numbers from 
Suntzeff et al. (1999). Our independent photometry of the 
comparison stars agrees well with their results. 

To measure the brightness of SN 1998bu, we performed 
aperture photometry of the supernova and comparison 
stars on each image. We derived an aperture correction 
determined from one or a few isolated bright stars mea- 
sured through circular apertures of varying radii (with sky 
measurements in a surrounding annulus). We were thus 
able to measure the supernova light through a small aper- 
ture, so that noise from the background sky would be min- 
imized. Due to the varying seeing and pixel scales at the 
sites, we did not impose a fixed angular size aperture for 
all the observations. In all cases, though, the aperture 
for the supernova instrumental magnitudes was the same 
size as the aperture for the comparison stars in the field. 
We chose aperture photometry over point-spread-function 
(PSF) fitting primarily out of convenience, but also be- 
cause in many instances the fields were too small or the 
exposures too short to reliably determine the PSF from 
nearby stars. 

We also determined the linear color transformation coef- 
ficients for each telescope/instrument/filter combination. 
At sites which had photometric conditions, these were the 
same as those used in determining comparison star mag- 
nitudes. For the other sites, color terms were determined 
from observations of either Landolt fields or other standard 
star fields but allowing for a varying zero point to account 
for the non-photometric conditions. The transformation 
coefficients were also checked with the local comparison 
stars. Color terms for each telescope/filter combination 
are listed in Table 0. Uncertainties in the transformation 
coefficients were propagated to the error estimate in the 
final photometry. 

Since our observations of SN 1998bu used local standard 
stars at the same airmass as the supernova, no airmass 
correction was necessary for the differential photometry. 
By using the measured color terms, we only had to deter- 
mine the zero point of each frame, by solving for the offset 
between the comparison star color-corrected instrumen- 
tal magnitudes and their standard magnitudes. The zero 
point was determined from the flux-weighted combination 
of all available comparison stars. The scatter in the com- 
parison star magnitudes (typically ^ 0.02 mag) was used 
as an estimate of the internal error in determining the 
zero point. From the derived zero point and color trans- 
formation coefficient, we were then able to transform the 
supernova instrumental magnitude into a standard mag- 
nitude. 

This procedure was quite effective in the majority of 
cases. However, supernova photometry requires special 
care. Where the color terms are large (due to a mismatch 
between the filter/detector response and a standard Lan- 
dolt response), the use of a linear color correction can be 
insufficient. This problem is especially acute in supernova 



photometry because of deep, wide features in the spectra 
of supernovae. Even when a large color correction works 
well when applied to stars (i.e., yields very small differ- 
ences with standard magnitudes), such large corrections 
may not be appropriate to supernovae. In most cases, fil- 
ters were well-matched and the derived color corrections 
were small, so nothing more complicated was required. 
However, in two data sets (the Gettysburg and Welles- 
ley B filters), the color-corrected supernova magnitudes 
were significantly discrepant with other data. In both of 
these cases the color terms were quite large. To combat 
this problem, we determined a correction based on the B 
filter transmission and detector sensitivity functions from 
these two sites. We used these response functions with 
spectrophotometric observations of a number of SNe la at 
varying ages, and determined magnitude corrections rela- 
tive to the standard passband defined by Bessell (1990). 
This procedure mirrors the use of K-corrections for high- 
redshift supernovae (Kim, Goobar, & Perlmutter 1996; 
Schmidt et al. 1998). Combining these corrections (usu- 
ally ~ 0.1 to 0.2 mag) with zero points determined from 
the comparison stars brought these two data sets into good 
accord with the others, though the photometry has signf- 
icantly higher uncertainty (typically ^0.1 mag). 

In addition, we encountered a puzzle in the MDM R- 
band data, which was discrepantly bright (by ~ 0.1 mag) 
compared to data from four other telescopes at nearly the 
same epoch. However, in this case, the derived R color 
term was not very large. We were unable to procure the fil- 
ter and detector characteristics for this data set. Thus we 
cannot provide an adequate explanation for this anomaly; 
it is possible that the filter transmission is mismatched 
just such that a color term derived from stellar observa- 
tions (which were somewhat limited in color) would be 
small (e.g., high transmission at some wavelengths com- 
pensated by low transmission at others), while the super- 
nova spectrum at the epoch of the observations led to a 
discrepant magnitude. Since four other telescopes pro- 
vided mutually consistent observations at the same epoch, 
we have disregarded the MDM R data. We were fortu- 
nate to possess contemporaneous observations so that we 
could discover this discrepancy, and it illustrates the need 
for a careful investigation of filter and detector character- 
istics when combining observations of supernovae made at 
various sites (Suntzeff et al. 1999; Wells et al. 1994). 

A futher complication in supernova photometry is 
proper discrimination between light from the supernova 
and light from the underlying galaxy (Boisseau & Wheeler 
1991). Accurate subtraction of the galaxy background is 
essential to measure correct magnitudes, and even more 
important in measuring light curves, as a constant un- 
removed galaxy background will cause a supernova light 
curve to look more shallow (i.e., mimic a slower decline 
rate). Though SN 1998bu is projected on a spiral arm of 
M96, the background from the galaxy is relatively faint. 
We were fortunate to possess UBVRI images of M96 prior 
to the appearance of SN 1998bu, taken with the FLWO 
1.2-m for another program. Our original plan was to use 
this image as a template and apply the galaxy subtraction 
techniques described by Schmidt et al. (1998; see also Fil- 
ippenko et al. 1986). This was very successful in a few 
cases, but it turned out that the template image quality 
was much poorer (due to seeing and the pixel scale) than 




Fig. 1. — SN 1998bu in M96 with local comparison stars. The field is 7'6 x 7'6; north is up and east is to the left. 
The image is a composite of several y-band observations of SN 1998bu taken with the F. L. Whipple Observatory 1.2-m 
telescope in May of 1998. 



Table 1 
Local standard star UBVRIJHK magnitudes. 



star 



V 



R 



H 



K 



1 (S6) 


13.521(0.027) 


13.594(0.012) 


13.068(0.009) 


12.774(0.012) 


12.475(0.011) 


12.02(0.05) 


2(S7) 


15.437(0.029) 


15.551(0.015) 


15.016(0.013) 


14.715(0.015) 


14.403(0.015) 




3 (SI) 


15.523(0.035) 


15.520(0.018) 


14.895(0.017) 


14.553(0.019) 


14.199(0.019) 




4(S2) 


16.501(0.040) 


16.495(0.027) 


15.789(0.024) 


15.384(0.026) 


14.971(0.025) 




5 (S8) 


16.778(0.039) 


16.280(0.025) 


15.441(0.024) 


14.987(0.027) 


14.573(0.025) 




6 (S12) 




19.210(0.042) 


17.730(0.031) 


16.805(0.040) 


15.583(0.033) 


14.36(0.05) 



11.85(0.03) 11.77(0.03) 



13.75(0.04) 13.54(0.03) 



Note. — The S identifiers are from SuntzefF et al. (1999). Uncertainties in the magnitudes are listed in parentheses. 



Table 2 
Photometric Color Terms 



Telescope 



U -B 



B-V 



V-R 



V ■ 



V{B - V) 



CfA4 

CfAl 

MDM 

KP09 

WIYN 

Gett 

Well 

KAIT 



0.935 
0.926 
0.782 
0.907 

1.140 



0.887 


0.973 


0.939 


0.982 


1.061 


0.98r 


0.912 


1.005 


0.974 


1.060 


0.549*^ 


0.894 


0.658'' 


0.980 


0.929 


1.190 



1.086 
1.068 
1.027 
1.027 
1.100 
1.006 
1.099 
1.082 



-0.041 
-0.027 
-0.019 
-0.009 
-0.019 
-0.009 
-0.030 
-0.001 



Note. — The telescope designations are: CfA4, FLWO 1.2-m with 4Shooter; CfAl, FLWO 1.2-m with AndyCam; 
MDM, MDM Hiltner 2.4-m; KP09, KPNO 0.9-m; WIYN, WIYN 3.5-m; Gett, Gettysburg 0.4-m; Well, Wellesley 0.6-m; 
and KAIT, KAIT 0.76-m. The tabulated values for the color columns {U-B,B~V,V-R, V-I) are the transfor- 
mation coeffiencents from the standard color to the instrumental color, e.g., CfA4: {u — h) ~ 0.935(C/ — i?) + const. The 
tabulated values for the last column are transformation coefficients from B — V to the instrumental v magnitude, e.g., 
CfA4: v-V = +0.041 (B -V)^ const. 

''The MDM i?-band photometry was not used; see text for details. 

''As described in the text, we have not used a linear transformation to place Gettysburg and Wellesley i3-band supernova 
photometry onto the standard system. 



Table 3 
UBVRI Photometry of SN 1998bu. 



Julian Day 


U 


B 


V 


R 


/ 


Telescope 


2450944.68 






12.45(0.14 


12.23(0.10) 


11.95(0.15 


KAIT 


2450947.63 




12.47(0.11) 


12.15(0.02 


11.86(0.04) 


11.72(0.04 


Gett 


2450948.59 






12.08(0.02 


11.80(0.04) 


11.68(0.04 


Well 


2450948.65 




12.41(0.09) 


12.09(0.02 


11.80(0.04) 


11.70(0.04 


Gett 


2450949.67 






11.99(0.03 


11.74(0.05) 


11.67(0.06 


Gett 


2450949.67 


11.96(0.04) 


12.29(0.02) 


11.98(0.01 


11.72(0.03) 


11.63(0.03 


CfA4 


2450951.58 




12.28(0.07) 


11.90(0.02 


11.71(0.04) 


11.68(0.04 


Well 


2450951.63 






11.93(0.03 


11.72(0.04) 


11.67(0.04 


Gett 


2450951.67 






11.91(0.03 


11.80(0.09) 




KAIT 


2450951.68 




12.21(0.02) 


11.90(0.01 


11.68(0.03) 


11.66(0.03 


CfA4 


2450952.62 




12.29(0.07) 


11.87(0.02 


11.70(0.04) 


11.70(0.04 


Well 


2450952.64 






11.88(0.02 


11.71(0.05) 


11.74(0.05 


Gett 


2450952.66 


12.00(0.05) 


12.21(0.02) 


11.86(0.01 


11.66(0.03) 


11.71(0.03 


CfA4 


2450952.67 






11.85(0.03 


11.75(0.09) 


11.69(0.06 


KAIT 


2450953.63 




12.34(0.11) 








Gett 


2450953.70 


12.03(0.05) 


12.24(0.02) 


11.87(0.01 


11.66(0.03) 


11.73(0.03 


CfA4 


2450955.64 


12.14(0.05) 


12.28(0.02) 


11.86(0.01 


11.65(0.03) 


11.79(0.03 


CfA4 


2450955.66 




12.26(0.11) 


11.88(0.02 


11.66(0.04) 


11.80(0.05 


Gett 


2450955.68 






11.85(0.03 




11.82(0.06 


KAIT 


2450956.59 




12.35(0.07) 


11.90(0.02 


11.71(0.04) 


11.83(0.04 


Well 


2450956.64 


12.23(0.04) 


12.32(0.02) 


11.87(0.01 


11.67(0.03) 


11.80(0.03 


CfA4 


2450956.68 






11.89(0.03 


11.79(0.09) 


11.84(0.06 


KAIT 


2450957.63 


12.31(0.05) 


12.36(0.02) 


11.92(0.01 


11.72(0.03) 


11.87(0.03 


CfA4 


2450957.65 




12.34(0.11) 


11.96(0.03 


11.75(0.05) 


11.88(0.05 


Gett 


2450959.68 


12.48(0.04) 


12.50(0.02) 


12.01(0.01 


11.82(0.03) 


11.98(0.03 


CfAl 


2450960.57 




12.57(0.07) 


12.02(0.02 


11.91(0.04) 


12.04(0.04 


Well 


2450960.69 


12.61(0.05) 


12.56(0.02) 


12.04(0.01 


11.89(0.03) 


12.05(0.03 


CfAl 


2450961.57 




12.74(0.11) 


12.12(0.03 


12.02(0.05) 




Gett 


2450961.60 




12.64(0.07) 


12.07(0.02 


11.99(0.04) 


12.10(0.04 


Well 


2450962.66 


12.76(0.05) 


12.70(0.02) 


12.12(0.01 


12.05(0.03) 


12.21(0.03 


CfA4 


2450963.64 


12.86(0.05) 


12.81(0.03) 


12.21(0.02 


12.15(0.03) 


12.26(0.04 


CfA4 


2450963.68 




12.88(0.05) 


12.26(0.02 


12.28(0.06) 


12.30(0.04 


KAIT 


2450964.63 


13.01(0.05) 


12.94(0.03) 


12.28(0.01 


12.21(0.03) 


12.30(0.03 


CfA4 


2450964.68 




12.97(0.05) 


12.34(0.02 




12.32(0.04 


KAIT 


2450965.64 


13.14(0.05) 


13.03(0.02) 


12.35(0.01 


12.29(0.03) 


12.31(0.03 


CfA4 


2450965.67 




13.05(0.03) 


12.42(0.02 




12.33(0.04 


WIYN 


2450966.60 




13.10(0.07) 


12.42(0.02 


12.35(0.04) 


12.29(0.04 


Well 


2450967.65 


13.45(0.05) 


13.24(0.02) 


12.48(0.01 






CfA4 


2450967.66 


13.35(0.05) 


13.24(0.02) 


12.50(0.01 




12.35(0.03 


MDM 


2450968.67 


13.53(0.04) 


13.40(0.02) 


12.59(0.01 


12.46(0.03) 


12.33(0.03 


KP09 


2450968.68 


13.50(0.04) 


13.36(0.02) 


12.55(0.01 




12.33(0.03 


MDM 


2450968.72 


13.60(0.05) 


13.43(0.03) 


12.62(0.02 


12.48(0.03) 


12.31(0.03 


WIYN 


2450969.59 




13.51(0.11) 


12.62(0.02 


12.43(0.05) 


12.26(0.04 


Gett 


2450969.60 




13.40(0.07) 


12.61(0.02 


12.43(0.04) 


12.25(0.04 


Well 


2450969.65 


13.68(0.04) 


13.51(0.02) 


12.64(0.01 


12.47(0.03) 


12.31(0.03 


KP09 


2450969.67 


13.75(0.05) 


13.52(0.03) 


12.66(0.02 


12.50(0.03) 


12.30(0.03 


WIYN 


2450969.69 




13.53(0.05) 


12.66(0.02 


12.50(0.06) 


12.32(0.06 


KAIT 


2450969.69 


13.66(0.04) 


13.47(0.02) 


12.59(0.01 




12.28(0.03 


MDM 


2450970.69 


13.82(0.05) 


13.59(0.02) 


12.64(0.01 




12.25(0.03 


MDM 


2450971.71 


13.97(0.05) 


13.72(0.02) 


12.69(0.01 




12.23(0.03 


MDM 


2450972.65 


14.12(0.04) 


13.84(0.02) 


12.74(0.01 




12.20(0.03 


MDM 


2450972.68 


14.16(0.04) 


13.89(0.02) 


12.80(0.01 


12.52(0.03) 


12.22(0.03 


WIYN 


2450973.58 




13.95(0.11) 


12.80(0.02 


12.45(0.05) 


12.16(0.04 


Gett 


2450973.68 




14.01(0.05) 




12.54(0.06) 


12.22(0.04 


KAIT 


2450974.67 


14.45(0.05) 


14.11(0.03) 


12.87(0.02 


12.51(0.04) 


12.19(0.04 


WIYN 


2450978.69 




14.51(0.06) 


13.09(0.02 


12.61(0.06) 


12.12(0.05 


KAIT 



many of our observations, and degrading the supernova 
observations to match the template added undesirable cor- 
related noise. In additional cases the observed fields did 
not align well with the template image (because there were 
very few stars in some of the small- field observations). In- 
stead we estimated the galaxy background flux directly 
from the template image (using an aperture and sky an- 
nulus appropriate to each observation) and subtracted this 
flux from that of the supernova. The correction to the su- 
pernova magnitude was initially negligibly small and grew 
larger as the supernova faded, but the maximum correc- 
tion made was only 0.025 mag. Our prediscovery images 
justify the assumption of a small host-galaxy flux contribu- 
tion made by Suntzeff ct al. (1999). Our flnal photometry 
for SN 1998bu is fisted in Table | and tfic optical light 
curves are shown in Figure 0. 

2.3. Infrared Photometry 

Infrared photometry of SN 1998bu in the JHK pass- 
bands was obtained at the Fred L. Whipple Observatory 
(FLWO), the Mount Stromlo Observatory (MSO), the In- 
frared Telescope Facility (IRTF), and witfi tfie ESO New 
Tecfinology Telescope (NTT). The FLWO data were ob- 
tained with the 1.2-m telescope and "STELIRCam" IR 
camera which consists of two 256^ InSb detector arrays 
permitting simultaneous imaging in two filters (ToUestrup 
et al. 1999). The FLWO filters were manufactured by 
Barr Associates in 1987 for a number of institutions in- 
cluding NOAO. The MSO data was taken with the 2.3-m 
telescope and CASPIR Spectrograph/Imager which uses a 
256^ InSb detector. The IRTF data was obtained with the 
NASA 3-m telescope and NSFCAM IR camera which uses 
a 256^ InSb detector. The NTT is a 3.5-m aperture tele- 
scope and observations were made using the SOFI imaging 
spectrograph which employs a 1024^ HgCdTe array. 

Flat fields and sky frames were created using offset field 
images staggered between the supernova exposures. The 
alternating offset field frames were subtracted from the 
corresponding data images and then divided by the nor- 
malized, average fiat field. Standard stars for the FLWO 
observations were taken from Elias et al. (1982), while 
the NTT used standards from Persson et al. (1998), and 
the IRTF used the UKIRT faint standards (Hunt et al. 
1998). The MSO data was calibrated with standards from 
Carter & Meadows (1995). Where possible, the results 
were transformed to the Elias system, but the variety of 
detectors, filters and standards introduces systematic er- 
rors on the order of 0.05 mag. Two stars (stars 1 and 
6) near the supernova were calibrated from the FLWO 
data and used as secondary standards on non-photometric 
nights; their magnitudes are listed in Table pi. The galaxy 
background at the position of the supernova is smooth 
and much fainter than the supernova light, so that aper- 
ture photometry was sufficient. The resulting supernova 
photometry is given in Table U and displayed in Figure ^ 

2.4. Optical Spectroscopy 

All of our optical spectroscopic observations were ob- 
tained at the FLWO 1.5-m telescope with the FAST spec- 
trograph. This long-slit spectrograph has been designed 
for high throughput and features a thinned, back-side il- 
luminated, anti-reflection coated CCD detector. The slit 
length is 180" and can be adjusted to several widths; we 



have generally employed a 3" slit for our observations. 
These were made using a 300 line/mm grating which re- 
sults in a resolution of roughly 0.6 nm and a usable wave- 
length range from 360 to 720 nm. We began observations 
of SN 1998bu on May 16.1 and continued through July 
15.1, on a total of 27 nights. 

We have reduced the spectra in the standard manner 
with IRAF. The two-dimensional CCD exposures were 
corrected for bias and dark-current and were flat-flelded 
using CCDPROC. A ID spectrum was extracted at the 
supernova position subtracting the neigboring sky using 
the APEXTRACT task. Wavelength calibration was per- 
formed by extracting the same aperture from an exposure 
of a HeNeAr lamp taken just after the supernova obser- 
vation and identifying emission lines. We also performed 
flux calibration through the reduction of a spectrophoto- 
metric standard star each night (Massey et al. 1988). The 
conditions for the supernova observations were usually not 
photometric, so the absolute flux calibration is generally 
unreliable. The relative flux measurement may also suffer 
because of differential refraction (Filippenko 1982), as the 
slit was always oriented east-west (PA = 90°) rather than 
being optimally oriented at the parallactic angle. In gen- 
eral, then, our spectra likely underestimate the fiux in the 
blue by 10 to 20%. In addition, since we have not used a 
blocking filter, second-order light contaminates the red end 
of the spectra. The supernova flux in the blue is generally 
lower than in the red and the detector sensitivity to the 
blue photons is also low, so that the second-order contam- 
ination is on the order of only a few percent. Second order 
contamination in the standard star spectra is more signifi- 
cant (since the standard stars are typically quite blue), so 
that the flux calibration in the red is also somewhat un- 
certain. Uncertainties in the flux calibration in both the 
blue and the red thus make these spectra unsuitable for 
spectrophotometry. Contamination of the supernova spec- 
tra by underlying galaxy light (after subtracting the local 
sky) was small, as in the case of the optical photometry. In 
some cases, multiple observations on the same night were 
combined into one. Cosmic rays and telluric lines were 
removed by hand. Figure shows a representative subset 
of our optical spectroscopy of SN 1998bu and the spectral 
evolution of the supernova covering 60 days, from approx- 
imately 3 days before maximum light in the B band. The 
complete spectroscopic data set is available upon request. 



2.5. Infrared Spectroscopy 

Infrared spectra of SN 1998bu were obtained with the 
KPNO 4-m telescope and OSU-NOAO Infrared Imager 
(ONIS) on 1998 June 14.2 UT. The ONIS covered the 
iiT-band range from 2.0 /im to 2.4 /im with a resolution 
of 1.6 nm (FWHM). A sequence of two-minute exposures 
was made, while stepping the target along the slit. These 
were combined for a total integration of 24 minutes. A 
spectrum of the F5V star BS 4281 was divided into the 
supernova spectrum to remove telluric absorption (except 
in the deep absorption bands at 1.4 and 1.9 /im). A smooth 
spectrum of an F5V star was created by interpolating the 
broad band colors (Johnson et al. 1968), with the zero 
point set by the catalog magnitude of BS 4281. Multiply- 
ing the supernova spectrum and the synthetic F5V star 
spectrum corrects for the detector sensitivity variations 



Table 3 — Continued 



Julian Day 


U 


B 


V 


R 


/ 


Telescope 


2450979.70 




14.60(0.06) 


13.14(0.03) 


12.64(0.06) 


12.14(0.05) 


KAIT 


2450980.65 


15.15(0.04) 


14.68(0.02) 


13.16(0.01) 


12.61(0.03) 


12.08(0.03) 


CfA4 


2450981.66 


15.22(0.04) 


14.74(0.02) 


13.25(0.01) 


12.69(0.03) 


12.12(0.03) 


CfA4 


2450981.69 




14.74(0.06) 


13.24(0.03) 


12.71(0.07) 


12.14(0.05) 


KAIT 


2450982.60 






13.30(0.03) 


12.79(0.06) 


12.13(0.05) 


Gett 


2450982.65 


15.29(0.05) 


14.81(0.02) 


13.30(0.02) 


12.75(0.03) 


12.16(0.03) 


CfA4 


2450982.70 




14.84(0.06) 


13.36(0.02) 


12.83(0.06) 


12.24(0.04) 


KAIT 


2450983.66 


15.31(0.05) 


14.89(0.02) 


13.38(0.01) 


12.83(0.03) 


12.24(0.03) 


CfA4 


2450984.67 


15.39(0.04) 


14.94(0.02) 


13.45(0.01) 


12.92(0.03) 


12.33(0.03) 


CfA4 


2450984.70 




14.95(0.05) 


13.51(0.02) 


13.00(0.06) 




KAIT 


2450985.65 




15.00(0.02) 


13.53(0.01) 


13.03(0.04) 


12.42(0.03) 


CfA4 


2450986.66 


15.48(0.05) 


15.05(0.02) 


13.56(0.01) 


13.08(0.03) 


12.47(0.03) 


CfA4 


2450986.69 




15.09(0.05) 




13.13(0.06) 




KAIT 


2450987.66 


15.50(0.04) 


15.09(0.02) 


13.63(0.01) 


13.12(0.03) 


12.55(0.03) 


CfA4 


2450988.66 


15.53(0.04) 


15.12(0.02) 


13.67(0.01) 


13.18(0.03) 


12.63(0.03) 


CfA4 


2450989.65 


15.58(0.04) 


15.18(0.02) 


13.72(0.01) 


13.24(0.03) 


12.68(0.03) 


CfA4 


2450990.65 


15.61(0.05) 


15.21(0.02) 


13.77(0.01) 


13.28(0.03) 


12.75(0.03) 


CfA4 


2450990.69 




15.22(0.05) 


13.82(0.02) 


13.34(0.06) 


12.80(0.05) 


KAIT 


2450991.65 


15.62(0.05) 


15.22(0.02) 


13.82(0.01) 


13.31(0.03) 


12.81(0.03) 


CfA4 


2450992.64 


15.64(0.04) 


15.24(0.02) 


13.84(0.01) 


13.38(0.03) 


12.86(0.03) 


CfA4 


2450993.64 


15.62(0.06) 


15.26(0.03) 


13.88(0.02) 


13.42(0.03) 


12.91(0.03) 


CfA4 


2450994.64 


15.69(0.05) 


15.29(0.02) 


13.92(0.01) 




12.97(0.03) 


CfA4 


2450995.64 


15.69(0.04) 


15.31(0.02) 


13.96(0.01) 


13.49(0.03) 


13.04(0.03) 


CfA4 


2450996.64 






13.97(0.01) 




13.09(0.03) 


CfA4 



Note. — The telescope designations are as in Table g. 



Table 4 
JHK Photometry of SN 1998bu. 



Julian Day 


J 


H 


K 


Telescope 


2450945.6 


11.76(0.06) 


11.88(0.06) 


11.81(0.05) 


FLWO 


2450948.6 


11.59(0.06) 


11.77(0.06) 


11.59(0.05) 


FLWO 


2450951.88 


11.66(0.04) 


11.84(0.04) 


11.60(0.03) 


IRTF 


2450970.7 


13.32(0.06) 


11.94(0.05) 


11.95(0.05) 


FLWO 


2450974.91 


13.23(0.06) 


11.68(0.06) 


11.89(0.05) 


MSO 


2450975.97 


13.12(0.04) 


11.79(0.03) 


11.77(0.03) 


ESO 


2450976.88 


13.08(0.06) 


11.65(0.06) 


11.77(0.05) 


MSO 


2450978.65 


12.81(0.06) 


11.73(0.06) 


11.74(0.05) 


FLWO 


2450978.87 




11.67(0.10) 


11.77(0.10) 


MSO 


2450984.76 


12.68(0.05) 


12.00(0.03) 


12.05(0.04) 


IRTF 
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Fig. 2.— UBVRI light curves of SN 1998bu. 
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Fig. 3.— JHK light curves of SN 1998bu. 
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Fig. 4. — Optical spectra of SN 1998bu labelled by epoch with respect to B maximum. For clarity the spectra have been 
shifted vertically by arbitrary amounts. Unrcmoved telluric features are marked. 
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and roughly calibrates supernova flux. The spectral flux 
was then adjusted to match the observed iiT-band pho- 
tometric magnitude of the supernova interpolated to the 
date of the spectrum. 

A spectrum was also taken with the SOFI instrument 
on the ESO NTT on 1998 June 11.0 UT. SOFI covered 
0.95 /xm to 2.5 /im in two grating settings with significant 
overlap. Four 120-s exposures were obtained at each tilt at 
four slit positions allowing good sky subtraction. Spectra 
of HD177619 (an F7V star) were used to remove telluric 
bands and calibrate the relative sensitivity of the detector 
as described above. Absolute flux calibration was done us- 
ing the _ff-band magnitude determined from SOFI imaging 
done on the same night. 

The KPNO and ESO spectra were combined into a sin- 
gle high-quality spectrum of SN 1998bu at an age of about 
-1-25 days which is shown in Figure H. There are few good 
infrared spectra of SNe la at the same epoch to compare 
with these data; however, a spectrum of SN 1995D taken 
with the MMT-FFSPEC at an age of +24 days is also 
shown for comparison. 

3. RESULTS 
3.1. Spectra 

Type I supernovae are distinguished from those of type 
II by the absence of hydrogen in their spectra, and SNe 
la are further distinguished from SNe lb and SNe Ic by 
the prominent Si II (A635.5 nm) absorption at maxi- 
mum light (for examples and a detailed discussion see 
Filippenko 1997). Other absorption features in the op- 
tical at maximum light are predominantly from interme- 
diate mass elements (Si, Ca, S, O, Mg) at high veloc- 
ity (~ 10,000 kms~^) in the outer layers of the super- 
nova ejecta. Lines of Fe become prominent at about two 
weeks after maximum light as the effective photosphere re- 
cedes into the ejecta, and about a month after maximum 
light the supernova enters the optically-thin nebular phase 
where forbidden emission lines of iron-peak elements (Fe, 
Co) dominate (Kirshner & Oke 1975). The optical spectra 
of SN 1998bu shown in Figure ^ follow this typical evolu- 
tion. A more detailed comparison is illustrated in Figure]^ 
where we show the optical spectra of SN 1998bu and other 
prototypical SNe la near maximum light. The spectra are 
remarkably similar, though there are some differences in 
the detailed shapes and velocities of the features. 

A more quantitative comparison between SN 1998bu 
and other typical SNe la is illustrated in Figure 0. We 
show the velocities of Si II A635.5 nm and Ca II H & K 
flux minima as a function of supernova phase for a number 
of prototypical SNe la: SN 1994D (Patat et al. 1996), SN 
1992A (Kirshner et al. 1993), SN 1990N (Leibundgut et 
al. 1991a), SN 1989B (Barbon et al. 1990; Wells et al. 
1994), and SN 1981B (Branch et al. 1983). SN 1998bu 
falls well within the scatter defined by these other objects. 

However, the measurement of Ca II H & K velocities is 
made more challenging by the presence of a feature blue- 
ward of the Ca feature (indicated in Figure ^ by a ques- 
tion mark), which may be due to either Si or possibly 
high-velocity Ca (Branch 1998, personal communication; 
Hatano et al. 1999). For SN 1998bu this unidentified fea- 
ture was well separated from the normal Ca feature, as 
is seen in the early spectra of the sequence in Figure m 



The feature weakens with time such that by day -1-28 only 
the normal Ca feature remains. The sequence suggests 
that the Ca absorption velocity should be defined by the 
redder of the two troughs. For other SNe la. Figure || 
shows that the unidentified feature and the Ca feature are 
not always well separated, such as in the SN 1989B and 
SN 1990N spectra, where only a shoulder is visible rather 
than two distinct minima. The spectrum of SN 1972E was 
taken at significantly lower resolution, and there it looks 
like a single absorption feature. Thus comparing measured 
Ca II H & K velocities is tricky; our use of the red trough 
may lead to a systematically lower velocity measured at 
early times, as seems to be the case in Figure M. 

Supernovae with unusual luminosity also have unusual 
spectra. Both SN 1991T, one of the most luminous SNe la 
known, and SN 1991bg, one of the least luminous, showed 
spectral peculiarities in their maximum light spectra. SN 
1991T lacked a well-defined Si II A635.5 nm feature at 
maximum (Filippenko et al. 1992a; Phillips et al. 1992), 
though its subsequent evolution was similar to other SNe 
la. On the other hand, maximum light spectra of SN 
1991bg showed a deep absorption trough attributed to 
Ti II centered at a wavelength of 420 nm (Filippenko et 
al. 1992b; Leibundgut et al. 1993). The absence of spec- 
troscopic pecuHarities in SN 1998bu makes it a suitable 
calibrator of the SN la distance scale. 

The infrared spectrum of SN 1998bu in Figure || is one 
of the best obtained for a type la supernova. It is strik- 
ingly similar to an IR spectrum of SN 1995D taken at the 
same phase. SN 1995D also showed no spectroscopic pe- 
culiarities in the optical and a typical light curve (Riess 
et al. 1999), bolstering the evidence that SN 1998bu is a 
fairly normal SN la. The IR spectra are qualitatively a 
good match to the model spectra of Wheeler et al. (1998), 
although the models do not extend to the observed age 
of SN 1998bu and are meant to fit the peculiar event 
SN 1986G. In the A'-band {^ 2.2 /im), which consists 
of absorption features of Co, Ni, and Si, the spectra of 
SNe la 1998bu, 1995D and 1986G (Wheeler et al. 1998; 
Frogel et al. 1987) are nearly identical from 14 days past 
maximum onwards. The iJ-band (~ 1.6 /im) spectra of 
SN 1998bu and SN 1995D do not show as large a gap be- 
tween the peaks at 1.6 /im and 1.8 /im as does SN 1986G. 
The deficit at 1.7 //m is not visible at all in the peculiar 
SN 1991T (Bowers et al. 1997), but that spectrum was 
taken at a more advanced age than the others. These ob- 
servations suggest that the 1.7 /im gap depth may possibly 
be correlated with light curve decline rate (and therefore 
luminosity) and supports the idea of Wheeler et al. (1998) 
that the 1.7 /im gap is an indicator of the highest velocity 
of the Ni/Co region, though clearly more infrared spectra 
of SNe la are required to test this hypothesis. As in the 
optical spectra, the infrared spectra do not indicate that 
SN 1998bu was peculiar. 

3.2. Light Curves and Peak Brightness 

The optical light curves of SN 1998bu presented in Fig- 
ure are among the best sampled early-time light curves 
of any SN la. The U band observations are particularly 
valuable as SNe la have not often been observed in this 
passband. In addition, our observations began before max- 
imum light (unless stated otherwise, we take maximum 
light to mean the time of maximum brightness in the B 
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band). These light curves are typical of SNe la; the UBV 
curves are well fit by the templates of Leibundgut (1989). 
In section 3.3 we describe the detailed analysis of these 
light curve shapes, an essential part of using this SN la as 
a distance indicator. 

Work on SNe la as standard or calibrated candles em- 
ploys the maximum brightness in the B and V bands. We 
have determined the peak apparent magnitude and time 
of maximum in the i?-band for SN 1998bu using a simple 
quadratic fit to the points within roughly five days of the 
light curve peak, weighted by their photometric uncertain- 
ties. The results are 



JDBmax = 2450952.8 ±0.8, 
^Bmax = 12.22 ± 0.03. 



(1) 
(2) 



(Throughout the paper, all apparent and absolute magni- 
tudes are expressed in units of mag). The supernova ap- 
parent magnitude in V at the time of B maximum given 
above is 

VBmax = 11.88 ±0.02. (3) 

The time when the supernova was brightest in the V 
band was JD 2450954.4 ± 1.0, at V = 11.86 ± 0.02. These 
results are quite consistent with the finding of Leibundgut 
(1989) that maximum light in V occurs about two days 
after maximum light in B and his result that the V mag- 
nitude at that time is 0.02 mag brighter than the V mag- 
nitude at the time of B maximum. The time of maximum 
brightness in the U band is not well determined because of 
the starting point of our data set, but our data are consis- 
tent with the Leibundgut (1989) result of U maximum oc- 
curring three days before maximum light in B. The U ap- 
parent magnitude at B maximum is t/Bmax = 12.01 ±0.05. 
SNe la show increased variety in their R and / light curves 
as compared to the bluer passbands, illustrated clearly in 
the composite light curves presented by Riess et al. (1999). 
The R and / light curves of SN 1998bu are consistent with 
these composite curves both in their general shape and 
times of maximum in those bands. 

In comparing our photometry with the independent data 
set of Suntzeff et al. (1999), we find excellent agreement. 
The times of maxima in the various passbands are quite 
consistent given the stated uncertainties, as are the light 
curves in general. In Figure g, we compare the photometry 
directly by plotting the difference between our magnitudes 
and those of Suntzeff et al. (1999). To make this compar- 
ison we have spline-interpolated the Suntzeff et al. (1999) 
light curves to the times of our observations, and have only 
computed differences when observations were within two 
days of each other to ensure the accuracy of the interpo- 
lation. The uncertainty in the magnitude difference was 
taken as the quadrature sum of the stated photometric 
uncertainties. As in the case of the photometry itself, the 
difference uncertainties are significantly correlated (due to 
the uncertainty in the comparison stars, for instance) . The 
largest differences occur, as expected, in the L'^-band, and 
in general the agreement between the two data sets is best 
near maximum light. Given the difficulties particular to 
supernova photometry, the consistency in the light curves 
is reassuring. Because small systematic differences in pho- 
tometry can have a magnified effect in distance determi- 
nation (through the reddening, for example), proper ac- 



counting of the (correlated) photometric errors is vital in 
order to obtain consistent results. 

One important result from the photometry is that 
the observed color of the supernova at maximum, {B — 
^)Bmax — 0.34 ± 0.04 mag, is significantly redder than 
typical SNe la, which have {B-V)Bma.x — 0.00 ±0.04 mag 
(Schaefer 1995). Very underluminous supernovae such as 
SN 1991bg can have quite red intrinsic colors at maxi- 
mum, but they also show distinct spectroscopic peculiari- 
ties. The absence of any such peculiarities in SN 1998bu 
suggests that this red color is not intrinsic to the super- 
nova but rather a result of interstellar extinction along 
the line of sight. Supporting evidence for this extinction 
is presented in Section ^. 

The optical color curves of SN 1998bu are also quite 
typical, as shown in Figure |[ where we present the color 
evolution of SN 1998bu compared to SN 1989B, also a 
spectroscopically normal SN la (Wells et al. 1994). Both 
of these supernovae had a similar B — V color index at 
maximum light and the general shape of the color evolu- 
tion of these two supernovae are in reasonable agreement. 
The slope of the B — V rise is measurably different, re- 
sulting from the fact that SN 1989B was a faster-declining 
object. The other striking difference is the offset of the 
U — B curves for the two supernovae. This could be a 
result of photometric uncertainties in calibrating the U 
band, where detector sensitivities and filter transmissions 
can differ substantially from one site to another and re- 
quire careful calibration (see the discussion by Suntzeff et 
al. 1999). However, it may also point to interesting diver- 
sity in the C/-band characteristics of SNe la, or diversity 
in the selective-to-total extinction properties of dust. A 
more detailed investigation of SNe la light curves in the 
near ultraviolet is warranted. This may be particularly im- 
portant for observations of SNe la at high redshift where 
observations at optical wavelengths probe the rest-frame 
ultraviolet. Without a thorough understanding of SNe la 
C/-band properties, cosmological inferences based on rest- 
frame [/-band light curves are suspect. 

In the infrared, the light curves of SN 1998bu match the 
JHK ternplates developed by Elias et al. (1985), as shown 
in Figure 0, where we have fit the templates to the data by 
adjusting them independently in magnitude and together 
in time. The bright second maximum typical of SNe la in 
the infrared passbands shown in the templates is clearly 
observed in SN 1998bu, as is the J-band deficit. The 
best-fit magnitude offsets to the templates are as follows: 
J = 12.14±0.09, H = 11.99±0.05, and K = 12.04±0.05. 
We find that the fiducial time to as defined by Elias et al. 
(1985) is about 3 days before maximum light in B. This 
differs slightly from the Elias et al. result, which suggested 
io was roughly 5 days before maximum light; however, 
only three supernovae were used in that determination, so 
it would not be surprising if there were some variance. A 
larger sample of infrared light curves, especially with ob- 
servations near the first maximum, would be useful. We 
note that Meikle & Hernandez (1999) show a large amount 
of infrared photometry of SN 1998bu near optical maxi- 
mum; combining these data with our light curve (which 
is dominated by later points) should yield one of the best 
infrared light curves of a type la supernova. 

3.3. Multicolor Light Curve Shape Analysis 
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The relation between luminosity and light curve shapes 
for SNe la, as quantified by Phillips (1993), and subse- 
quently Hamuy et al. (1996a), led to the development 
of techniques to measure distances to SNe la from multi- 
color light curves. Hamuy et al. (1996b) showed how BVI 
light curves and templates (Hamuy et al. 1996d) could be 
combined to derive accurate distances using a x^ analysis. 
In a similar vein, Riess, Press, & Kirshner (1996a) devel- 
oped the Multicolor Light Curve Shape (MLCS) method, 
a statistical technique to measure distances to SNe la from 
their BVRI light curves, allowing for reddening in the host 
galaxy. In this approach, the light curves of a "training 
set" of supernovae with estimated luminosities and extinc- 
tions are used to derive template light curves for a fiducial 
SN la, along with derived correction templates which de- 
tail the change in the light curve shapes as a function 
of luminosity and extinction. We focus on MLCS dis- 
tances to SNe la in this paper; Phillips et al. (1999) 
present an extension to their template-fitting technique 
which also incorporates reddening, and the results of ap- 
plying this method to SN 1998bu are reported by Suntzeff 
et al. (1999). 

The orginal MLCS training set was based on nearby SNe 
la and relative distances measured to their host galaxies 
via the Tully-Fisher (TF), the surface brightness fluctua- 
tion (SBF), or the planetary nebula luminosity function 
(PNLF) methods. The only requirement was that these 
methods give accurate relative distances to the galaxies. 
Once trained, the method can be used on the light curves 
of a SN la, to determine the luminosity difference. A, be- 
tween that supernova and the fiducial (A = 0) supernova 
as well as a derived extinction to the supernova. Appli- 
cation of MLCS to a sample of more distant supernovae 
indicated the effectiveness of this procedure. The disper- 
sion in SN la distances about the Hubble line was reduced 
from a ~ 0.4 mag in the standard candle assumption to 
a ~ 0.12 mag with MLCS. The technique also demon- 
strated the Hubble law was applicable to large distances 
corresponding to velocities cz ~ 30,000 kms"""^, as was 
also shown by Hamuy et al. (1996b). 

However, uncertainties in these secondary distance de- 
terminations, inherent difficulties in estimating the extinc- 
tion to some supernovae, and the small number of training 
set objects caused some problems in the derived MLCS SN 
la distance scale, as pointed out by Saha et al. (1997). To 
combat the major problems, Riess et al. (1998a) presented 
a new version of MLCS in which the relative distances for 
the training set objects were derived from their host galaxy 
recession velocities and the Hubble law for galaxies with 
redshifts cz > 2500 km s~^. In addition, estimates of the 
extinction to the training set supernovae were refined, and 
effects of extinction on the shape of the light curve based 
on temporal variations in the effective selective-to-total ex- 
tinction ratios from the evolving supernova spectrum (Nu- 
gent, Kim, & Perlmutter 1999) were also included. The 
procedure detailed in Riess et al. (1998a) was restricted to 
B and V light curves up to 40 days past maximum, to be 
applied to the high-redshift supernovae that are the focus 
of that paper, but a procedurally identical version (with 
one exception) using full BVRI light curves is used in this 
paper. Here we have used a "wide" choice of the a priori 
extinction distribution, with cr(yly) ~ 1 mag rather than a 
distribution which overly simplified the models of Hatano, 



Branch, & Deaton (1998). We found that our distribu- 
tion based on those models predicts many fewer signifi- 
cantly reddened supernovae than are observed. Applying 
a Bayesian filter based on that particular distribution leads 
to underestimated extinctions in reddened supernovae and 
produces biased distances. We have reverted to the less re- 
strictive prior distribution that was used by Riess, Press, 
& Kirshner (1996a), so that the posterior distribution is 
shaped primarily by the observations rather than by the 
prior (cf . Figure |l2|) . 

The MLCS analysis fits the light curves with templates 
that are a function of luminosity and extinction. To deter- 
mine the peak brightness of the supernova (in the V band 
which is arbitrarily chosen as the default), rather than us- 
ing just the points near maximum light, the whole light 
curve is used, through a weighted average of the difference 
between the light curve and the best-fit template. We 
designate this weighted average (Vemax) to differentiate it 
from the estimate of maximum light brightness based on 
only the points near the time of maximum light, Vsmax- 
Again, we note that the time of maximum light is defined 
in the B band, such that both Vemax and (Vemax) describe 
the V magnitude at the time of peak B luminosity. The 
difference between Vemax and (Vemax) is quite small in all 
cases, but disregarding the distinction would make it ap- 
pear as though there are discrepancies in derived distances 
at the level of ~ 0.02 mag. 

The MLCS analysis of a set of supernova light curves 
results in three major parameters: (Vemax), ^y, and A. 
Apphcation of the MLCS method on the BVRI light 
curves of SN 1998bu yields a maximum light brightness 
(Vemax) = 11.89, which is quite consistent with the result 
from just the points near maximum, Vemax = 11.88±0.02. 
The derived extinction is Ay = 0.94 mag and the lumi- 
nosity difference A = 0.02 mag (i.e., the supernova was 
2% less luminous than the fiducial). For the purposes of 
comparing supernovae and measuring distances, we define 
my = (Vemax) — ^y as the extinction-corrected maximum 
light apparent magnitude of the supernova in the V band. 
We further define the quantity rrty = (Vemax) — Ay — A, 
which would be the maximum light apparent V-band 
brightness of the supernova had it been free of absorp- 
tion and of fiducial luminosity. For SN 1998bu, then, we 
have 



my = 10.95 ±0.18 and 



m?, = 10.93 ±0.18, 



(4) 

(5) 



where the uncertainties are primarily due to the (corre- 
lated) uncertainties in the derived luminosity correction 
and extinction, with the uncertainty in the extinction 
(±0.15 mag) being the dominant component. 

The derived extinction can be compared directly with 
the supernova's red color at maximum light. Assuming 
that a typical unreddened SN la has {B — V)emax — 
0.00 ± 0.04 (Schaefer 1995) implies a color excess for SN 
1998buof S(S-V) = 0.34±0.06mag. Adopting (at max- 
imum light) Rv = 3.1 yields Ay = 1.05 ± 0.19 mag, fully 
consistent with the MLCS extinction derived from the full 
BVRI light curves. We note that the expected extinc- 
tion from our Galaxy along the line of sight to SN 1998bu 
is small, E{B — V) = 0.025 mag (Schlegel, Finkbeiner, 
& Davis 1998), so that the bulk of the reddening is from 
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M96 itself! The B - V color of SNe la is also gener- 
ally quite uniform at ~ 35 days past maximum light, with 
B - V ~ 1.1 ± 0.1 mag (Lira 1995; Riess et al. 1998a; 
Phillips et al. 1999). The observed color of SN 1998bu at 
that time, B — V = 1.48 ± 0.04 mag, also implies a color 
excess consistent with the extinction derived from the full 
MLCS analysis. Suntzeff et al. (1999) derive a total red- 
dening for SN 1998bu of E{B - V) = 0.37 mag based on 
the B - y and V-I color evolution (Phillips et al. 1999), 
which is consistent with our results. 

Further evidence comes from the infrared light curves, 
where effects of dust are expected to be small. By compar- 
ing our IR observations with those of other well-observed 
type la events and assuming that the optical-IR colors are 
constant for normal SN la, we can derive another check on 
the inferred extinction. For example, the V maximum of 
SN 1981B was fainter by 0.08 mag than the V maximum 
of SN 1998bu, while in iiT-band the SN 1981B light curve 
was fainter than the SN 1998bu light curve by 0.87 mag 
(Elias et al. 1981). Using the extinction law of Cardelli, 
Clayton, & Mathis (1989) with Ry = 3.1, the difference in 
visual extinction between the two supernovae is given by 
AAv = (A7^- Ay)/0.886. Hence, the relative extinction 
between 1981B and 1998bu is AAy = 0.89 mag. Unfortu- 
nately, there are few type la supernovae with low extinc- 
tions and good IR light curves which can be used in this 
way. Combining infrared and optical data of SN 1980N 
and SN 1981D, both in NGC 1316 (Hamuy et al. 1991), 
gives a relative extinction of 0.81 mag and a rough lower 
limit consistent with our derived value. SN 1989B was 
highly extinguished as is SN 1998bu. The visual magni- 
tude difference between these two is 0.11 mag (SN 1998bu 
is brighter) and the K difference is —0.07 mag (SN 1998bu 
is fainter), suggesting that there is 0.20 mag less visual ex- 
tinction to SN 1998bu than SN 1989B. Wells et al. (1994) 
found a color excess to SN 1989B of E{B -V) = 0.37 mag, 
meaning the extinction to SN 1998bu would be Ay = 0.95 
mag. All of the estimates from the infrared photometry 
are consistent with a total visual extinction to SN 1998bu 
of Ay = 0.9 ±0.2 mag, supporting the value derived in the 
MLCS analysis. 

These checks based on the color excesses at a number of 
wavelengths from B-band to iiT-band are consistent with 
the view that the shape of the extinction curve is likely 
close to the standard Galactic law (Riess, Press, & Kirsh- 
ner 1996b), though the absolute normalization is not con- 
strained. Spectrophotometric observations of SNe la have 
been used to determine the normaHzation, with Ry ~ 3.1 
at maximum light, as well as the temporal variation of the 
effective Ry due to the evolution of the supernova spec- 
trum (Nugent, Kim, & Perlmutter 1999). The t/-band 
photometry still provides some cause for concern; if the 
blue U — B color is due to differences in the properties 
of the absorbing dust, the estimated extinction may be 
incorrect. 

Independent checks of the MLCS-derived extinction are 
valuable. One such check is the presence of interstellar 
Na I D and Ca II H & K absorption in high-resolution 
spectra, which was reported by Munari et al. (1998) and 
Centurion et al. (1998). The equivalent width of the Na I 



Dl (A589.0 nm) absorption lines reported by Munari et 
al. (1998) were 0.019 nm and 0.035 nm at velocities cor- 
responding to our Galaxy and M96, respectively, and we 
confirm these measurements even with our low-dispersion 
spectra. Using the Munari & Zwitter (1997) calibration of 
the correlation between the equivalent width and redden- 
ing they derive color excesses of 0.06 and 0.15 mag, for a 
total reddening of E{B — V) = 0.21 mag. However, the re- 
lation between the absorption-line equivalent widths and 
the reddening has a large scatter, with a typical disper- 
sion of 0.15 mag in E{B — V) for multi-component lines 
(Munari & Zwitter 1997), so that these values do not con- 
tradict the extinction inferred from the MLCS analysis (cf. 
discussion by Suntzeff et al. 1999). 

We can also check our derived luminosity with other 
techniques. The relation between light curve shape and 
luminosity was pioneered by Phillips (1993) and Hamuy 
et al. (1995, 1996a) using the quantity Ami5(i3), which 
parameterizes the BVI light curves in terms of the B mag- 
nitude decline of the supernova over the fifteen days af- 
ter maximum light. From our light curve, we measure 
Ami5(i3) = 1.02 ± 0.04, which agrees very with Suntzeff 
et al. (1999), who found Arrii^iB) = 1.01 ± 0.05. Di- 
rect comparison of the luminosity correction is made dif- 
ficult because of effects of extinction on the light curve 
shape, as well as differences in our respective fiducial tem- 
plates. Nevertheless, we can measure Ami5{B) from the 
MLCS fiducial template (A ^ 0, Ay = 0) which yields 
ATOi5(i3) = 1.08. Thus the measured values of Amis and 
A, both of which imply SN 1998bu to be quite close to the 
MLCS fiducial template, indicate consistency in the two 
approaches. This is not surprising, since both methods use 
the shape of the observed light curve in a similar fashion. 

An alternative approach was described by Nugent et al. 
(1995), who presented correlations between spectral fea- 
tures and intrinsic SNe la luminosity. In particular they 
define two indicators: 7?.(Si II), the ratio of the depths of 
Si II absorption lines observed at 580 and 615 nm; and 
TZ{Ca II), the flux ratio of the continuum levels just blue 
and red of the Ca II H & K absorption. Our maximum- 
light spectra of SN 1998bu yield 7^(Si II) = 0.23±0.02 and 
7^(Ca II) = 1.47 ± 0.04. Riess et al. (1998b) present Hn- 
ear relations between A and both TZ{Si II) and TZ{Ca II), 
which yield a mean luminosity correction for SN 1998bu of 
A = -0.02±0.16 and A = 0.06±0.22, respectively. These 
agree well with the luminosity correction result derived in 
the MLCS analysis. 

All indications thus suggest that SN 1998bu was an in- 
trinsically normal type la supernova, significantly extin- 
guished by dust along the line of sight. With the MLCS 
analysis we determine the peak brightness ((Vemax) = 
11.89), the luminosity correction (A = 0.02 mag), and 
the extinction {Ay = 0.94 mag), so that SN 1998bu can 
be used to calibrate the SN la distance scale. 

4. THE DISTANCE SCALE 

4.1. Huhhle-Flow SNe la 

Constructing the Hubble diagram requires a sample of 
well-observed SNc la in the Hubble flow where errors due 

^Stanok (1998) and Arce & Goodman (1999) have recently concluded that the Galactic reddening maps of Schlegel, Finkbeiner, & Davis 
(1998) might overestimate the extinction in regions where E{B — V)> 0.15 mag. The Galactic extinction towards M96 is well below this level, 
so this should not be a major concern. 
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to peculiar velocities are expected to be small, and which 
is analyzed in exactly the same way as the local calibra- 
tors. Our MLCS sample consists of 42 SNe la, 26 from 
the Calan/Tololo supernova search (Hamuy et al. 1993, 
1996c) and 16 from the CfA supernova monitoring cam- 
paign (Riess et al. 1999). The only further selection crite- 
ria we have imposed (other than those inherent in the two 
data sets) is a cut in the host-galaxy recession velocity, 
which has been corrected to the cosmic microwave back- 
ground (CMB) frameQ. We have excluded supernovae in 
galaxies with cz < 2500 km s~^ where peculiar motions 
become increasingly important. Additionally, we have ex- 
cluded supernovae in galaxies with logcz [km s""'^] > 4.5, 
where the relation between luminosity-distance and red- 
shift begins to be non-linear at a level which could affect 
our results (Schmidt et al. 1998). 

The Hubble-flow sample consists of the following SNe 
la: SN 1990O, SN 1990af, SN 1991U, SN 1991ag, SN 
1992J, SN 1992K, SN 1992P, SN 1992ae, SN 1992ag, SN 
1992al, SN 1992aq, SN 1992au, SN 1992bc, SN 1992bg, SN 
1992bh, SN 1992bk, SN 1992bl, SN 1992bo, SN 1992bp, 
SN 1992br, SN 1992bs, SN 1993B, SN 1993H, SN 19930, 
SN 1993ac, SN 1993ae, SN 1993ag, SN 1993ah, SN 1994M, 
SN 1994Q, SN 1994S, SN 1994T, SN 1995D, SN 1995E, 
SN 1995ac, SN 1995ak, SN 1995bd, SN 1996C, SN 1996Z, 
SN 1996bl, SN 1996bo, and SN 1996bv. We emphasize 
that all of these supernovae have well-sampled multicolor 
CCD light curves, with photometry obtained, reduced, and 
transformed to the standard system in a similar fashion, 
an essential feature which allows us to combine the data 
sets. 

We use MLCS to turn these SNe la into standard can- 
dles, so that the extinction- and luminosity-corrected max- 
imum light absolute magnitude. My = my — /z, is a con- 
stant, where n is the distance modulus. Using the defini- 
tion of the distance modulus, 

ml- - M^- = ^j = 5 log d + 25 = 5 log ^ + 25, (6) 

Ho 

where d is the distance in Mpc, cz is the recession velocity 
in km s^^, and Hq is the Hubble constant measured in its 
conventional units of km s~^ Mpc""'^, we get the result 



logcz — 0.2my = logiJo 



0.2M^ 



(7) 



Since the absolute magnitude of the fiducial SN la is taken 
to be constant, we can determine that 



logcz — 0.2™" 



ay. 



(8) 



Here ay = logHo — 0.2AIy — 5 is "the intercept of the ridge 
line" and is a constant which can be determined from ob- 
servations of Hubble-flow SNe la alone. 

In Figure HO we present this relation for our sample 
of 42 Hubble-flow SNe. The small scatter allows us to 
solve precisely for the intercept, ay, as shown more clearly 
in Figure O. The best-fit intercept (which is simply 



the mean of log cz 



0.2m'^j 



ay = 0.6772 ± 0.0049, where the uncertainty is the formal 
standard error in the mean, and assumes the residuals from 
each supernova are normally distributed and independent. 
We assign a ±300 km s""'^ la uncertainty to the redshift 
to account for the contribution of peculiar motions (Riess, 
Press, & Kirshner 1996a). The dispersion about the mean 
is (7{aY) — 0.0317; this corresponds to a dispersion in mag- 
nitudes (obtained simply by multiplying by 5) of ct = 0.16 
mag, implying only an 8% relative distance uncertainty per 
objectn. Our derived intercept depends upon the choice 
for the fiducial luminosity (A = 0) supernova. As long as 
the comparison with the local calibrators is made with the 
same choice, there is no problem. However, when com- 
paring this set of Hubble-flow SNe with those analyzed 
by a different technique, such as the Ami^{B) method 
(Suntzeff et al. 1999; Phillips et al. 1999), there will likely 
be an offset in ay due to the different choices of a fidu- 
cial luminosity. Such an offset will also be rcfiected as the 
same offset for the fiducial absolute magnitude, My. 

The quoted statistical uncertainty in the measurement 
of ay is quite small, and it is surely underestimated. 
We assumed that each supernova distance is independent, 
whereas in reality there exists some covariance. Thus, the 
true uncertainty in the mean does not simply decrease as 
a/y/N, but rather levels off due to a floor caused by sys- 
tematic uncertainties. It is thus important to estimate at 
what level this floor is reached. The formal uncertainty 
in ay corresponds to ±0.025 mag or just over one percent 
uncertainty in the Hubble constant arising just from the 
Hubble-flow supernovae. 

Sample differences between the Hubble-flow and cali- 
brating SNe la are one potential source of systematic un- 
certainty at this level. The present sample is necessarily 
imperfect; one difference occurs due to the host-galaxy 
type: our Hubble-flow sample includes supernovae in both 
early-type and late-type hosts, while hosts of the local cal- 
ibrators have Cepheid distances, and thus are of late-type 
only. Hamuy et al. (1995) have shown that SNe la in E/SO 
galaxies are systematically fainter than those in spirals or 
irregulars (we refer to galaxies classified Sa or later, in- 
cluding irregulars, as "spirals" in what follows). More ac- 
curately, the highest luminosity SNe la are found only in 
spirals (Riess et al. 1999), perhaps implying a relation be- 
tween recent star formation and the brightest SNe la. Nev- 
ertheless, SNe la brightness in both early-type and late- 
type galaxies correlate similarly with light curve shape, so 
that an MLCS (or Am^i^^B)) analysis will still correct SNe 
la in spirals and ellipticals to the same fiducial luminosity 
without using any information on the galaxy morphology, 
as demonstrated by Schmidt et al. (1998). To further test 
this we have divided our sample into two subsets, early- 
type and late-type hosts, solving independently for the 
intercept. The results are presented in Table ||. We see 
that after application of MLCS the offset between ay in 
early-type and late-type galaxies is inconsistent with zero 
at only the 1.5ct level. Determining whether this differ- 
ence is significant will require a larger sample. However, 



using all the supernovae is 

^Heliocentric redshifts for the host galaxies were first transformed to the Local Group rest frame by adding (-30, 297, -27) km s~^ in Galactic 
Cartesian coordinates (de Vaucouleurs et al. 1991; Lynden-Bell & Lahav 1988). The recession velocities in the Local Group frame were then 
transformed to the CMB rest frame by adding (10, -542, 300) km s"! (Smoot et al. 1992). 

^Some of this uncertainty arises from the uncertainty in the redshift due to peculiar velocities. Our sample has an error- weighted mean 
redshift of roughly 10,000 km s~^ so that the adopted 300 kms~^ uncertainty corresponds to a 3% distance uncertainty on average. This 
means the actual relative distance uncertainty intrinsic to the MLCS analysis of the supernovae is only 7% per object. 
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Fig. 10. — Apparent magnitude-redshift relation for 42 Hubble-flow SNe la which were corrected for extinction 
and to a fiducial luminosity with a multicolor light curve shape (MLCS) analysis. The supernovae are from the 
Calan/Tololo (Hamuy et al. 1996c) and CfA (Riess et al. 1999) data sets. Supernovae in late-type (>Sa) galaxies 
are shown with circles, those in early-type (E/SO) galaxies are shown with triangles. The best-fit ridge line is shown, 
logcz — 0.2my + 0.6772(±0.0049). The dispersion about the best-fit line is cr = 0.16 mag. 
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Fig. 11. — Ridge-line intercept versus recession velocity for the Hubble-fiow sample (left) and ridge- line intercept histogram 
{right). The best-fit mean intercept and its formal uncertainty are ay — 0.6772 ± 0.0049. 
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Table 5 
Intercepts of the Ridge Line for Hubble-Flow SNe Ia. 



Sample 



N 



ay 



a [mag] 



All 



42 



0.6772 ±0.0049 



0.16 



Late- type (>Sa) 
Early-type (E/SO) 

cz > 7000kms"i 
cz < 7000kms-i 



25 
17 

28 
14 



0.6716 ±0.0067 
0.6854 ±0.0067 

0.6712 ±0.0059 
0.6892 ±0.0081 



0.17 
0.14 

0.16 
0.15 



this offset would lead to a difference in the derived Hub- 
ble constant (for a fixed My) of ^ 3% (in the sense that 
early- type galaxies yield the slightly higher value). 

Another potential source of systematic uncertainty in 
the measurement of ay is the effect of galaxy peculiar 
velocities and flows. We have transformed the measured 
host-galaxy recession velocity to the frame at rest with re- 
spect to the CMB frame, but distortions of the velocity 
field will result in errors in our derived intercept. Some 
pecuHar velocity studies (e.g., GiovaneUi et al. 1998) 
have indicated convergence of the fiow field relative to 
the Local Group with the CMB dipole at redshifts cz ~ 
4000 km s~^, though others do not (e.g., Lauer & Postman 
1994). Our full sample is cut at cz > 2500 km s~^, so flows 
may be important at the lowest redshifts. To check this we 
examined a subset of data with cz > 7000 km s~^ where 
such motions should have a much smaller effect. As Table 
a shows, there is no signifcant difference in ay between our 
full sample and the sample restricted to cz > 7000 km s~^. 

There is a larger difference in comparing the more dis- 
tant sample with the remaining SNe la, i.e. those with 
cz < 7000 kms^^. In this case the offset in ay differs 
from zero at 1.8cr. Zehavi et al. (1998) have interpreted 
this result as the effect of a local void, whereby we live in a 
slightly underdense region compared to the average density 
of the universe, leading to a larger nearby expansion rate. 
For a fixed My the difference in ay would lead to a "lo- 
cal" Hubble constant larger than the global value by ~ 4%. 
Again, more Hubble-flow SN la light curves will help de- 
termine whether this offset is real, a statistical fluke, or an 
artifact of the analysis technique. Since the large majority 
of our sample (28 of 42) have cz > 7000 km s^^, the mean 
is more reflective of the global value. 

We have performed an additional test of our derived in- 
tercept by employing a simple geometric flow model which 
includes the effects of nearby mass concentrations such as 
the Virgo Cluster, the Great Attractor, and the Shapley 
Supercluster on the velocities of the supernovae host galax- 
ies (Mould et al. 1999). Using the position and redshifts 
of our galaxy sample, this model predicts that we underes- 
timate the Hubble constant by ~ 2% in assuming that the 
Hubble-flow galaxies are at rest with respect to the CMB 
frame. 

Additional sources of uncertainty may remain, for in- 
stance due to correlations in the MLCS analysis arising 
from the training set and construction of the templates, or 



uncertainties in the calibration of the photometric system. 
A larger sample of Hubble-flow SNe la would be helpful to 
determine what unidentified systematics may remain and 
at what level they affect our conclusions. 

Given that the identified systematic errors in the 
Hubble-fiow SNe do not yet definitively suggest a bias 
in our derived intercept, we use these results only as a 
guide to the size of the systematic uncertainty. Based on 
these explorations, a reasonable estimate of the la sys- 
tematic uncertainty in the Hubble-flow supernovae would 
be ±3% in the Hubble constant. Our best estimate for 
the intercept of the ridge-line and its total uncertainty 
is then ay — 0.6772 ± 0.0120. As we discuss below, 
this uncertainty in the Hubble-flow SNe la is dwarfed 
by both statistical and systematic uncertainties in the 
Cepheid-calibrated supernovae and our measurement of 
M^ (Hamuy et al. 1996b; Riess, Press, & Kirshner 1996a). 

4.2. Cepheid- Calibrated SNe la 

To measure the Hubble constant we need both the inter- 
cept of the ridge line, ay, and the maximum light absolute 
magnitude of our fiducial SN la. My . We use a sample of 
SNe la in galaxies with distances measured via Cepheids. 
SN 1998bu in NGC 3368 is only the most recent example, 
others are listed by Saha et al. (1997) in their Table 6: 
SN 1895B and SN 1972E in NGC 5253, SN 1937C in IC 
4182, SN 1960F in NGC 4496A, SN 1981B in NGC 4536, 
and SN 1990N in NGC 4639. Their table also includes 
SN 1989B in NGC 3627, but the distance to this galaxy is 
only inferred from distances to other galaxies in the Leo 
group (including NGC 3368). We restrict our sample to 
the best cases: supernovae in galaxies whose distances are 
directly measured by Cepheids rather than distances to 
groups or clusters. Schaefer (1998) has also recalibrated 
the fight curve of SN 1974G in NGC 4414, for which the 
host-galaxy Cepheid distance has been measured by the 
HST Distance Scale Key Project (Turner et al. 1998). 

The supernovae we use as calibrators are only those 
which are measured and analyzed in the same way as our 
Hubble-fiow sample, to avoid systematic errors. For this 
reason, we only consider SNe la which have photoelectri- 
cally measured multicolor light curves. As our observa- 
tions of SN 1998bu have shown, even CCD data measured 
and calibrated similarly can yield discrepant photometry 
depending on the details of the telescope, detector, and 
filters. Such calibration problems, as well as problems of 
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galaxy background subtraction and transformation to the 
Landolt system, make photographic photometry of super- 
novae subject to systematic differences in the peak bright- 
ness, colors, and the light curve shape. The last point is 
crucial; we must be certain that the observed light curve 
shape is an intrinsic property of the supernova, and not an 
artifact of systematic errors resulting from photographic 
photometry (Boisseau & Wheeler 1991; Pierce & Jacoby 
1995). While heroic efforts (e.g., Schaefer 1998 and refer- 
ences therein) have been made in compiling and reanalyz- 
ing older photographic (and even visual) light curves, the 
best path to the Hubble constant lies along another route: 
precise distances to well-observed objects. The drawback 
to a high standard for the data is that our calibrating 
sample is small, consisting of four SNe la: SN 1998bu, SN 
1990N, SN 1981B and SN 1972E. We pay in random error 
what we avoid in systematic bias and we believe this to be 
a good bargain. 

The Cepheid distances to the host galaxies of these four 
supernovae have been measured by two HST programs, 
but the general approach among all the HST Cepheid pro- 
grams is the same. From the derived mean magnitudes 
of the Cepheids and a PL relation (Madore & Freedman 
1991), distance moduh (relative to the LMC) can be deter- 
mined. The two-color photometry allows for an estimate 
of the extinction, either on a Cepheid-by-Cepheid basis 
or in the mean, yielding an extinction-corrected distance 
modulus for the host galaxy. 

It is important to use distance moduli for the host galax- 
ies which are derived consistently (e.g., always using the 
same PL relation, the same LMC distance modulus, the 
same extinction prescription, etc.) so that each supernova 
is on an equal footing. Thus we have not simply used 
the final distance modulus quoted in the Cepheid papers, 
but rather have tried to extract as uniform a set of dis- 
tance moduli as we can with limited information. This 
approach also allows us to more easily consider systematic 
effects in the Cepheid distances. For instance, in this sec- 
tion we use distance moduli with the LMC distance fixed 
at /iLMC = 18.50 mag (Hereafter, all distance moduli will 
have implied units of magnitudes). We do not yet include 
the uncertainty in this value because that uncertainty is 
implicit in each host galaxy distance, and moreover it is 
perfectly correlated, such that the derived mean absolute 
magnitudes will suffer the same uncertainty. We postpone 
discussion and quantification of such systematic ("exter- 
nal" ) uncertainties to the next section. 

4.2.1. SN 1998hu in M96 (NGC 3368) 

For SN 1998bu, Tanvir et al. (1995) discovered seven 
Cepheids with well-determined light curves in M96. They 
derived an extinction-corrected distance modulus of /i = 
30.32 ± 0.12 with the uncertainty coming from the pho- 
tometric errors and the uncertainty (in the mean) of the 
fit to the PL relation. However, their photometry was 
not corrected for the WFPC2 "long/short" exposure ef- 
fect (Hill et al. 1998) which leads to V and / magni- 
tudes systematically too bright by 0.05 mag in long expo- 
sures such as those for Cepheid programs. We have cor- 
rected the distance modulus of M96 for this effect, yield- 
ing fi = 30.37 ± 0.12, with the quoted uncertainty being 
only the "internal" error. Combined with the extinction- 
corrected maximum light apparent magnitude from Equa- 



tion 0, we derive an extinction-corrected maximum light 
absolute magnitude for SN 1998bu, My = -19.42 ± 0.22. 
To determine the absolute magnitude of our fiducial SN 
la, My, we also include the derived luminosity difference 
A = 0.02 mag for SN 1998bu, which leads to M^ = 
-19.44 ±0.22. 

Our absolute magnitude for SN 1998bu is fainter than 
that reported by Suntzeff et al. (1999) who found My = 
— 19.63 ± 0.19. This offset arises from different estimates 
of the extinction. We have employed Ay = 0.94 mag de- 
rived from the MLCS analysis of the BVRI light curves, 
whereas Suntzeff et al. (1999) derive a total E{B — V) = 
0.37 mag, yielding Ay — 1.15 mag. This difference is, for 
better or worse, within the uncertainties, and may arise 
partly from an offset in the intrinsic colors in the light 
curve fitting methods, but some more careful comparisons 
between the two methods may be necessary. Nevertheless, 
it turns out that the derived Hubble constant is not very 
different in the two methods. 

4.2.2. SN 1990N m NGC 4639 

For SN 1990N we have performed an MLCS analysis on 
BVRI light curves from Lira et al. (1998), with the results 
presented in Table g. We have used the extinction cor- 
rected Cepheid distance modulus to NGC 4639 determined 
by Saha et al. (1997), fi == 32.03 ±0.22, based on 15 high- 
quality Cepheids. From this we derive an absolute magni- 
tude for SN 1990N oi My = -19.78±0.33 and an estimate 
of the fiducial absolute magnitude My = — 19.46 ± 0.33. 

4.2.3. SN 1981B in NGC 4536 

Our MLCS analysis of SN 1981B was based on the BVR 
hght curves of Buta & Turner (1983). Saha et al. (1996) 
found a total of 73 Cepheids in NGC 4536 and determined 
an extinction corrected distance modulus, fi = 31.10±0.13. 
This leads to an absolute magnitude for SN 1981B of 
My = -19.46 ±0.23 andM{^ = -19.12 ± 0.23. 

4.2.4. SN 1972E in NGC 5253 

Finally, for SN 1972E, we have BVI light curves 
from Ardeberg & de Groot (1973) and Leibundgut et 
al. (1991b). Saha et al. (1995) have presented their fi- 
nal analysis of Cepheids in NGC 5253, with their result 
^ = 28.08 ± 0.10. However, some cautions are in order. 
These observations were made with the original WFPC in- 
strument (with spherical aberration) rather than WFPC2; 
also, the J-band light curves were transformed from the 
HST F785LP filter rather than F814W. These differences 
may possibly lead to a small systematic difference in the 
derived distance modulus compared to the other host 
galaxies. Furthermore, the Cepheid sample is small, as 
there are only 5 Cepheids with high-confidence mean mag- 
nitudes in both V and /. Using the derived apparent 
moduli in V and /, and estimating the Cepheid extinc- 
tion in the same manner as for the other three host galax- 
ies, we have rederived the same distance modulus as Saha 
et al. (1995) but with a significantly larger uncertainty, 
/i = 28.08 ± 0.26. Our MLCS analysis of SN 1972E then 
leads to My = -19.80 ± 0.29 and M^ = -19.42 ± 0.29. 

We present the MLCS results, Cepheid distances, and 
absolute magnitudes for SN 1998bu and the three other 
calibrating SNe la in Table o. We note that the es- 
timates for My are consistent given their uncertainties, 
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though SN 1981B seems to give a measurably fainter value. 
Since the estimates are mutually consistent and there is 
no a priori reason to distrust any of them, we take the 
data at face value. The error-weighted mean gives the 
maximum light absolute magnitude of our fiducial SN la, 
My — — 19.34±0.13. As with the ridge-line intercept, here 
again we must be wary of the uncertainty estimate since it 
ignores the covariance in the Cepheid distances. Kochanek 
(1997) has shown there is significant statistical covariance 
in the Cepheid distance moduli, even beyond the common 
zero-point set by the LMC distance, which arises from a 
number of sources including the Cepheid photometry, the 
determination of mean magnitudes, and the fit to the PL 
relation. If we include an estimate of this statistical covari- 
ance to determine how to combine the Cepheid distances 
(usually incorrectly assumed to be independent) by using 
the Pearson correlation cocfRcient, r ~ 0.5 (C. Kochanek, 
personal communication), our best estimate of the fiducial 
absolute magnitude and its statistical uncertainty becomes 
My — — 19.34±0.17. This estimate still does not incorpo- 
rate some sources of systematic uncertainty, includi ng t he 
LMC distance, which we dicuss in detail in Section |5.2|. 



4.3. The Hubble Constant 

With estimates of the ridge-line intercept and the 
maximum-light absolute magnitude of our fiducial SN la, 
we derive the Hubble constant, from 



logiJo 



0.2M?/ 



ay- 



(9) 



The mean of the four calibrating SNe la gives My — 
-19.34 ± 0.17, and using aU 42 Hubble-flow SNe la gives 
av — 0.6772 ± 0.0120, which results in our best estimate 
of the Hubble constant. 



Hq = 64.4I^;5'kms"' Mpc" 



(10) 



where the uncertainty does not include systematic uncer- 
tainties in the Cepheid distance scale to be discussed be- 
low. Even so, it is important to note that the total uncer- 
tainty in ay (±0.0120) is about three times smaller than 
the statistical uncertainty in 0.2M{^(±0.034). The sta- 
tistical error in this small sample of calibrating SNe la, 
arising from both the uncertainty in the luminosity- and 
extinction-corrected supernova brightness and the uncer- 
tainty in the Cepheid distance moduli, dominates the sta- 
tistical uncertainty in Hq (Hamuy et al. 1996b; Riess, 
Press, & Kirshner 1996a). Reducing this statistical error 
can best be accomplished by observations of additional 
nearby SNe la and Cepheids in their host galaxies. SN 
1998bu is the first example where a new supernova has 
been studied in a galaxy where the Cepheid work is al- 
ready in the literature, but more will follow in the years 
ahead. The systematic uncertainty in the calibrating SNe 
and their Cepheid distances is still an important consider- 
ation, discussed in the next section. 

5. DISCUSSION 

5.1. Comparison with other work 

Suntzeff et al. (1999) used SN 1998bu and four other 
local calibrators (they included SN 1937C in IC 4182 in 
their local sample), combined with distant SNe la from 
the Calan/Tololo sample to derive a Hubble constant, 



Hq ~ 64.0 ± 2.2 km s^^ Mpc^^ (this uncertainty ignores 
covariances in the Cepheid distances). Their result, using 
a different method to convert light curves to luminosities 
and extinction (Phillips et al. 1999), is quite consistent 
with ours. There arc slight differences in measurements of 
individual objects, but the overall agreement is reassuring 
and indicates that SNe la are excellent distance indicators 
whose intrinsic diversity can be understood and quantified. 

The calibration of the peak absolute magnitude of SNe 
la has been driven by the great efforts of the HST pro- 
gram to measure Cepheid distances to the supernova host 
galaxies (Sandage et al. 1992). That group's latest pub- 
lished determination of the Hubble constant and statis- 
tical uncertainty (also ignoring Cepheid covariances) is 
iJo == 58±3 kms-^Mpc-i (Saha et al. 1997). This 
is consistent with our result ~ agreement to ~ 10% is 
good given the long history of measurements of the Hub- 
ble constant. Nonetheless, it is instructive to pinpoint 
where the differences arise. The Saha et al. (1997) analy- 
sis uses a "fiducial sample" of 56 Hubble-flow SNe la with 
B and V peak magnitudes, typically determined from pho- 
tographic plates. Enforced upon the sample is a velocity 
constraint, 3 < logcz[kms~^] < 4.5, and a color con- 
straint —0.25 < -Bmax — Knax < 0.20 mag (after correction 
for Galactic extinction) to avoid peculiar SNe la and those 
with large amounts of extinction. As local calibrators, they 
use seven SNe la in six galaxies with five Cepheid distances 
(the distance to SN 1989B was estimated by association 
to other Leo Group galaxies with measured Cepheid dis- 
tances; see their Table 6). Their local calibrator sample is 
not selected by the same criteria as their fiducial sample, 
as no color constraint was applied to the local calibrators. 
So SN 1895B was used with an estimated -Bmax but no 
color information, and SN 1989B was used though it was 
too red, with Bmax — ^ax = 0.35 ± 0.07 mag (Wells et 
al. 1994). (SN 1998bu would also be too red to meet the 
color requirement of their fiducial sample.) Nevertheless, 
excluding SN 1895B and 1989B from their analysis would 
only have a small effect and would increase their derived 
Hq by about 1 km s~^ Mpc^^. 

Though we have not used SN 1895B and SN 1989B in 
our local sample for reasons adduced earlier, there are ad- 
ditional differences between the two analyses. Saha et al. 
(1997) include no correction for variation in intrinsic SN la 
luminosity based on light curve shape. This is particularly 
significant for the calibrator sample where 3 out of the 4 
objects are slow decliners; including this correction makes 
the estimated fiducial absolute peak magnitude slightly 
fainter and explains about half our disagreement in H^. 
Most of the balance of the difference likely arises from the 
treatment of extinction. While Saha et al. (1997) correct 
some of their local calibrator peak magnitudes for extinc- 
tion individually, they do not apply an extinction correc- 
tion to their fiducial sample SNe. If the mean color excess 
of their fiducial sample were as little as E{B — ]/) ~ 0.03 
mag, correcting for extinction would increase the mean 
my by ~ 0.1 mag, and raise ay by ~ 0.02, which is the 
other half of the difference in Hq. Saha et al. (1997) argue 
that selection effects against the discovery of extinguished 
distant supernovae preclude a significant amount of extinc- 
tion in the fiducial sample. While this could be true, there 
is no demonstration that it is true for the sample they 
use. The selection effects in the several searches that led 



Table 6 
Cepheid-Calibrated SNe Ia 



27 



SN la 


Galaxy 


(VBmax) 


Av 


A 


O-MLCS 


/^Ccphcid 


Mv 


AfO 


1998bu 
1990N 
1981B 
1972E 


NGC 3368 
NGC 4639 
NGC 4536 
NGC 5253 


11.89 
12.68 
11.99 

8.43 


0.94 
0.43 
0.35 
0.15 


+0.02 
-0.32 
-0.34 
-0.38 


±0.18 
±0.25 
±0.18 
±0.13 


30.37 ±0.12 
32.03 ±0.22 
31.10 ±0.13 
28.08 ±0.26 


-19.42 ±0.22 
-19.78 ±0.33 
-19.46 ±0.23 
-19.80 ±0.29 


-19.44 ±0.22 
-19.46 ±0.33 
-19.12 ±0.23 
-19.42 ±0.29 


Mean 
















-19.34 ±0.17 



While it is important to understand the effects of metal- 
licity on Cepheid luminosities and colors in many wave- 
bands (see, e.g., Freedman & Madore 1990; Gould 1994; 
Stiff 1995), our particular concern is the effect on the HST 
Cepheid distances. Recently great efforts have been made 
to determine empirically the relation between metallicity 
and distance moduli measured in V and / with the stan- 
dard procedure of extinction correction. The "metallicity" 
dependence of extragalactic Cepheids is usually parame- 
terized in terms of [0/H], the logarithmic number abun- 
dance of oxygen to hydrogen, relative to solar composition 
and measured via H 11 region spectra. We adopt the no- 
tation 
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to the SNe la of their fiducial sample are quite complicated 
(see, e.g. Hamuy & Pinto 1999), and a priori statements 
about the possible extinction distribution of the distant 
supernovae are not, by themselves, evidence. In partic- 
ular, discovery of SNe la with Ay ~ 0.1 mag does not 
seem to be strongly suppressed. In Figure |2| we show the 
extinction distribution for our Hubble-flow and calibrat- 
ing samples and it is clear that a some of the supernovae 
found this way are, in fact, signifcantly extinguished. 

The MLCS analysis was designed to address these con- 
cerns; the supernova peak magnitude, extinction and lu- 
minosity correction are quantitatively estimated for each 
object, with a careful attention to correlations and the fi- 
nal distance modulus uncertainty. This obviates the need 
to make arbitrary sample cuts. We have taken particular 
care to analyze the local calibrators and the Hubble-flow 
SNe la by the same methods. While these technical differ- 
ences are significant, the difference in the derived Hubble 
constant between our approach and others is small. Uncer- 
tainty in the true Hubble constant still arises principally 
from the small size of the calibrator sample and from un- 
certainties in the lower rungs of the distance ladder. 

5.2. External uncertainties 

We examine three sources of systematic uncertainty in 
the distances of the galaxies that host the calibrating SNe 
la: the effect of metallicity on HST Cepheid distance mod- 
uli, recalibration of the Cepheid PL relation, and the dis- 
tance to the LMC, which are likely to be the most impor- 
tant sources of systematic error. 

5.2.1. Metallicity dependence of the Cepheid scale 

The remarkably tight PL relation of Cepheids in the 
LMC is quite a boon to distance measurement, but also 
calls for understanding variation in the PL (and therefore 
derived distance moduli) with environment. Metallicity, in 
particular, may play an important role. Theoretical stud- 
ies of Cepheid pulsation (e.g., Chiosi, Wood, & Capitanio 
1993) indicate metallicity can have an effect on the bright- 
ness and colors of Cepheids, though the size of this effect 
is uncertain. If the brightness and colors of Cepheids vary 
with metallicity, their distances will be misestimated, both 
because of an incorrect estimate of their intrinsic bright- 
ness as well as an incorrect estimate of the extinction based 
on the color excess. Systematic errors in distances would 
occur in Cepheid populations with significantly different 
mean abundances than the LMC. 

®The value of 'yyi = —0.44 quoted in their papers is for abundances measured in terms of [Fe/H]; 7v/ = —0.48 is correct for [O/H] (D. 
Sasselov, personal communication). 



which gives the change in the distance modulus per fac- 
tor of ten in metallicity, measured in mag/dex. Then the 
true distance modulus of a galaxy is given by /itruc = 
l^vi — 7y/([0/H] — [0/H]lmc), where ^vi is the mea- 
sured extinction-corrected distance modulus. Recently, 
Kennicutt et al. (1998) used HST observations of three 
fields in MIDI over which a large abundance gradient has 
been observed from measurements of H II regions. Their 
analysis showed a slight correlation of distance modulus 
with metallicity, with ^vi — —0.24 ±0.16. Beaulieu et al. 
(1997) and Sasselov et al. (1997) analyzed the effects of 
metallicity on the PL relation from observations of LMC 
and SMC Cepheids from the EROS microlensing project. 
Their conclusion was for a slightly stronger dependenceQ 
7v'/ = — 0.481;q'2. a global analysis of many Galactic and 
extragalactic Cepheids by Kochanek (1997) yielded a sim- 
ilar metallicity dependence; we consider this analysis in 
more detail in the next section. 

To assess the possible impact of metal abundance on 
the Hubble Constant from SNe la, we have recalculated 
the distance moduli to the calibrating galaxies with these 
empirically determined metallicity corrections. We adopt 
[O/H] - [0/H]lmc measurements for NGC 3368 and NGC 
5253 as compiled by Kennicutt et al. (1998) in their Ta- 
ble 4. The abundances for NGC 4639 and NGC 4536 
are from Kochanek (1997), who estimated these based 
on metallicity-magnitude and metallicity-galaxy type re- 
lations; we have adopted an uncertainty of ±0.20 on these 
estimates, though the results are not particularly sensi- 
tive to this choice. The effect of the recalculated distance 
moduli is shown in Table ul where we have derived My 
estimates for three cases: 7y/ = 0, i.e., no metallicity 
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Fig. 12. — Distribution of the total visual extinction, Ay, along the lines-of-sight to the 42 Hubble-flow SNe la (clear) 
and the 4 calibrators (shaded). The extinction estimates are derived from the MLCS analysis of the BVRI light curves. 
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dependence, 7y/ = —0.24 ± 0.16 from Kennicutt et al. 
(1998), and 7y/ = -0.48l°;^ from Beaulieu et al. (1997) 
and Sasselov et al. (1997). 

The results are interesting; the corrections for mctallic- 
ity yield little change in the inferred Hubble constant. Low 
Hq values from SNe la are not due to metallicity effects on 
the Cepheid distances to this sample of galaxies. This is 
primarily due to our inclusion of SN 1998bu in NGC 3368, 
which is the only metal-rich calibrator, by a factor of 0.7 
dex compared to the LMC. If metallicity were the culprit 
causing a low SNe la Hq, we would expect our inferred 
My for SN 1998bu to be very faint without a metallicity 
correction. However, the distance modulus of Tanvir et al. 
(1995) and our observations do not require a faint value. 
More strikingly, comparing SN 1998bu and SN 1972E we 
see that the estimates of My are almost exactly the same 
without a metallicity correction, even though the metal 
abundance is higher in NGC 3368 than in NGC 5253 by 
more than a factor 10. 

As a result, including a metallicity correction increases 
the dispersion in the estimates of My from the four cali- 
brating SNe la, from a{My) = 0.16 mag with no metal- 
licity correction, to a (My) = 0.21 mag for the Kennicutt 
et al. (1998) value, to (t(M^) = 0.29 mag for the Sasselov 
et al. (1997) value. The sample size of four is probably 
too small to place much confidence in this result, but if it 
were borne out by a larger sample, it would have interest- 
ing implications. Assuming that the intrinsic dispersion 
in the MLCS-corrected luminosity for the calibrating su- 
pernovae is the same as that for the Hubble-flow SNe la 
[a ~ 0.16 mag), any increase in the dispersion would arise 
from the Cepheid distance moduli. A significant increase 
in dispersion with a metallicity correction would imply ei- 
ther that the association between H II region metallicity 
and Cepheid metallicity is not straightforward, that the 
metallicity correction was incorrect, or that some other 
systematic uncertainty (perhaps in the properties of SNe 
in regions of different metallicity) was colluding with the 
metallicity to counter its effect in the uncorrected distance 
moduli. Increasing the sample of SNe la in galaxies with 
Cepheid distances (particularly covering a wide range of 
metallicity) would be a very desirable path to understand- 
ing this important uncertainty in the distance scale. We 
note that Nevalainen & Roos (1998) used this idea of "sta- 
tistical consistency" between Cepheid-calibrated distance 
indicators (including SNe la) to derive a metallicity depen- 
dence which brought the estimates of Hq into the best con- 
cordance, 7v/ — — 0.3ltoi4- Because of the high metal- 
hcity of NGC 3368, including SN 1998bu in their analysis 
would likely affect this result. 

In principle the supernovae themselves can provide an 
estimate of the metallicity correction, but minimizing the 
dispersion in our four calibrating supernovae yields jyi — 
-1-0.11 ± 0.37, which is only a weak constraint given the 
very small sample size. In addition, such a procedure ig- 
nores a possibly significant metallicity dependence in the 
brightness of the supernovae themselves. Thus the current 
supernova data do not provide strong evidence either for 
or against the incorporation of a metallicity dependence 
in the Cepheid distance moduli. We take the results at 
face value, and combined with the current best estimates 
of 7y/ we conclude that the systematic la error in Hq 
from metallicity considerations for this sample is small. 
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Mpc"^ 
5.2.2. The Cepheid PL relation 



We (and most authors) have adopted the F-band and I- 
band PL relations derived from LMC Cepheids by Madore 
& Freedman (1991). The relations are consistent with ear- 
lier estimates (e.g., Sandage & Tammann 1968; Feast & 
Walker 1987), but are based on a relatively small sample 
of objects compared to the number now known in more 
distant galaxies. An error in the PL slope or zero-point 
could be an important source of systematic uncertainty. 
In this section, by the zero-point we do not mean to in- 
clude the uncertainty in the distance to the LMC (still 
adopted as fihMC = 18.50, and discussed extensively in 
the next section), but only uncertainty in the zero-point 
that arises from a small sample of LMC Cepheids, even 
after assuming the LMC distance is perfectly known. 

The Madore & Freedman (1991) PL relations have sta- 
tistical zero-point uncertainties of 0.05 and 0.03 mag in 
V and / respectively. Unfortunately, residuals for the fit 
to the PL relations are correlated, and the PL slope is 
also slightly correlated to the zero point, so that the real 
uncertainty in the application of the PL relations is not 
just a straightforward quadrature combination of these un- 
certainties. Tanvir (1996) analyzed an augmented sample 
of LMC Cepheids with particular attention to correlated 
residuals and found a slightly tighter relation. However, 
his result indicated bias in the Madore & Freedman (1991) 
calibration such that for typical period distributions, HST 
Cepheid distance moduli are overestimated by ~ 0.1 mag, 
which would imply an increase in our estimate of My (to 
fainter intrinsic luminosity) by the same amount and a 
~ 5% increase in Hq for any given LMC distance. 

The global analysis of Kochanek (1997) carefully and 
consistently treats the Cepheid data to allow for distance 
estimates including the systematic uncertainties we have 
discussed so far. Again, the distance to the LMC is fixed 
at fi = 18.5. In Table H we present our estimates of My 
given the distance moduli derived by Kochanek (1997) and 
presented in his Table 3. We have added 0.05 mag to the 
NGC 3368 distance modului presented there to correct for 
the WFPC2 "long/short" exposure effect, which was not 
included in his analysis of the Cepheids in that galaxy (this 
is just a first approximation; the correct method would be 
to repeat his analysis with the fainter NGC 3368 Cepheid 
magnitudes). With that caveat, we have considered two of 
his models, the first being Model which derives distances 
in the "standard" method, with a global solution for the 
PL relation and correct treatment of correlated errors. As 
Table @ shows, the My estimates are generally fainter in 
this model, in line with the suggestion of Tanvir (1996), 
leading to a modest increase in Hq . The second model we 
consider is Model 3-15, which also derives a global solution, 
but further allows for effects of metallicity and positive ex- 
tinction. The result is again generally to decrease the host 
galaxy distances, and lead to a fainter My and higher Hq. 
However, we note that the estimates of My are not consis- 
tent with their uncertainties (which are derived from the 
quadrature sum of the uncertainty in my from MLCS and 
the uncertainty in the Cepheid distance modulus). Fur- 
thermore, the dispersion is much larger than the ~ 0.16 
mag expected from the SNe la alone; this implies that ei- 
ther the local calibrating SNe la are very different from 
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Table 7 
Effects of Metallicity 





Galaxy 


[O/H] - [0/H]lmc 




A/0 




SNIa 


■yvi = 


-yvi = -0.24 ±0.16 


7v/ = -0.481;;;^ 


1998bu 


NGC 3368 


+0.70 ±0.20 


-19.44 ±0.22 


-19.61 ±0.25 


-l9.78+_lfr 


1990N 


NOG 4639 


+0.10 ±0.20 


-19.46 ±0.33 


-19.48 ±0.33 


-19.5ltgf, 


1981B 


NGG 4536 


0.00 ±0.20 


-19.12 ±0.23 


-19.12 ±0.23 


-19.12+g:^l 


1972E 


NGC 5253 


-0.35 ±0.15 


-19.42 ±0.29 


-19.34 ±0.30 


-19.25l«:^o 


Mean 






-19.34 ±0.17 


-19.36 ±0.17 


-19.40 ±0.18 


Ho 






UAtli 


63.8+^-^ 


62.7l^:^ 



the Hubble-flow SNe la (a possibility we feel is unlikely 
based on their spectra and light curves), or that system- 
atic errors remain in the Cepheid distances derived from 
this particular model. Further analysis is required. 

These results suggest a la systematic uncertainty in the 
calibrating SNe la host galaxy Cepheid distance moduli 
of t^io '^^g due to the combined effects of metallicity, 
extinction and the calibration of the PL relation. (The 
resulting uncertainty in My is in the opposite sense, ^0^5 
mag.) Since definitive results regarding the exact mag- 
nitude of these effects are lacking, we just include these 



effects in our uncertainty, with My = —19.34 
estimate of the Hubble constant then becomes 



+0.20 

-0.18- 



Ho 



64:Atli km s" 



Mpc" 



Our 



(12) 



5.2.3. Distance to the LMC 



The final source of systematic error in the Hubble con- 
stant that we consider is the distance to the LMC. All 
methods to measure Hq which are based on HST Cepheid 
distances share this systematic uncertainty, so that com- 
parisons in the resulting Hq values between these meth- 
ods should not include an error component from the LMC 
distance (i.e., comparisons of HST Cepheid-calibrated H^ 

if 
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measurements should be to the analog of equation 
two such measurements disagree, they will disagree regard- 
less of the LMC distance). Formally, our best estimate of 
the Hubble constant, including the systematic uncertain- 
ties discussed above, is given by 



logiJo = 1.809 



+0.042 
0.038 



0.2 (AiLMC - 18.50) . 



(13) 



Of course, if we wish to compare our value of H^ with 
those derived from techniques which are independent of 
the LMC distance (e.g., SNe II expanding photospheres, 
gravitational lens time delays, Sunyaev-Zel'dovich effect, 
etc.), we have to provide a best estimate for hlmc and per- 
haps more importantly, its uncertainty. Measurement of 
quantities derived from the "true" Hubble constant, such 
as the dynamical age of the Universe, also requires such 
an estimate. 

However, measurements of the LMC distance modulus 
are notoriously inconsistent, ranging for the most part 
from - 18.2 to ~ 18.7. The value of /xlmc = 18.5 
that we have adopted has recently faced a strong chal- 
lenge from a "short" LMC distance based on Hipparcos- 
calibrated red clump stars (Stanek, Zaritsky, & Harris 



1998; Udalski 1998) and the study of detached eclipsing 
binaries such as HV2274 (Guinan et al. 1998), which 
give /iLMC = 18.18 ± 0.06 and 18.30 ± 0.07, respectively. 
However, recent applications of other methods, including 
Cepheids, RR Lyrae stars, Mira variables, and the SN 
1987A ring still yield a wide range of distance moduli, 
many with inconsistent error bars. As an exercise, we 
compiled a representative (though not exhaustive) sample 
of 19 LMC distance moduli published in the last two years 
(Alcock et al. 1997; Di Benedetto 1997; Feast & Catch- 
pole 1997; Gratton et al. 1997; Panagia et al. 1997; van 
Leeuwen et al. 1997; Whitelock, van Leeuwen, & Feast 
1997; Bergeat, Knapik, & Rutily 1998; Cole 1998; Fernley 
et al. 1998; Gieren, Fouque, & Gomez 1998; Gould & Uza 
1998; Guinan et al. 1998; Luri et al. 1998a; Luri et al. 
1998b; Madore & Freedman 1998; Oudmaijer, Groenewe- 
gen, & Schrijver 1998; Reid 1998; Udalski 1998). Naively 
assuming each measurement to be independent, a Bayesian 
analysis of these distance estimates in the spirit of Press 
(1997) yields a narrow probability density function (PDF) 
for the mean, with /ilmc = 18.55;|^Qg4, as shown in Figure 
[l3[ If we modify the analysis so that each distance method 
only gets one "vote" to reduce correlated errors, the PDF 
becomes quite asymmetric with /ulmc = 18.50to'i5i where 
the stated value is the peak of the PDF and the upper and 
lower uncertainties are derived from the points at which 
the cumulative probability is 0.841 and 0.159. This anal- 
ysis does not provide a reason for the discrepant values, 
and is subject to additional correlated errors that may still 
be lurking. Nevertheless, given the incompatible data, the 
method provides a reasonable and statistically defensible 
way to estimate ^lmCj and its uncertainty. 

As the current best estimate of the LMC distance mod- 
ulus, we adopt /j,lmc = 18.50lg']^5 niag. Others will un- 
doubtedly have differing estimates, and can use equation 



13 to determine the resulting Hubble constant. With our 
choice, we have as the estimate of the fiducial absolute 
magnitude, including all identified systematic uncertain- 
ties, M^ = -19.34l5]:2i- Thus, our final estimate of the 
Hubble constant incorporating this total systematic un- 
certainty is 



Ho = MAtli km s" 



Mpc" 



(14) 



5.3. Implications 



One direct implication of our derived Hubble constant 
is an estimate of the dynamical age of the universe, to, as- 



Table 8 
Results Based on the Global Analysis of Kochanek (1997) 
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SN la 



Galaxy 



Published'" 



M_l_ 



Model 0'^ 



Model 3-15'^ 



NGC 3368 


NGC 4639 


NGC 4536 


NGC 5253 



-19.44 ±0.22 
-19.46 ±0.33 
-19.12 ±0.23 
-19.42 ±0.29 



-19.40 ±0.33 
-19.54 ±0.41 
-18.99 ±0.28 
-19.04 ±0.35 



-19.54±0.23 
-19.21 ±0.28 
-18.63 ±0.22 
-18.74±0.20 



Mean 
Ho 



-19.34±0.17 

64.4t|j 



-19.20 ±0.22 



-19.03 ±0.27'= 

74 3+10-1 



'^Published refers to Aly estimates based on published distance moduli as described and modified in section 4.2 . 

''These two columns arc from Kochanek (1997). See text for details regarding the models. 

'^Since the values in this column seem to be inconsistent with their derived uncertainties, we have calculated an un- 
weighted mean. 
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Fig. 13. — Probability density functions for the mean LMC distance modulus based on a Bayesian analysis of 19 recent 
measurements. The dotted curve shows the pdf if each measurement is treated independently, while the solid curve shows 
the PDF when each distance measuring technique is given equal weight. For clarity, both distributions have been rescaled 
to peak at unity. 
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suming a Friedmann- Robertson- Walker cosmology. In an 
il]^j — 1, Einstein-de Sitter universe, HqIq — 2/3. How- 
ever, evidence from SNe la at high redshifts provides a 
strong observational constraint on this product, with Riess 
et al. (1998a) deriving Hoto = 0.93 ± 0.06 and Perlmutter 
et al. (1999) obtaining Hoto = 0.93 ±0.05. Assuming that 
systematic errors in these estimates are small (discussed 
extensively in both papers), we have 



^0 = 14.1 ± 1.6 Gyr. 



(15) 



This estimate is very nearly the same as that presented 
by Riess et al. (1998), and is in good accord with mea- 
surements of the ages of the oldest objects in the Universe 
(see Chaboyer 1998 for a review). The solution to the 
cosmological "age crisis" is not to be found in alternative 
estimates of Ho, but rather in discarding the ^Im = 1 cos- 
mology (which would require ip = 10.1 ±1.1 Gyr). 

The supernovae can also tell us about the structure of 
groups and clusters. For instance, Graham et al. (1997) 
obtained a Cepheid distance to NGC 3351 which is also 
a member of the Leo I group along with NGC 3368. 
However, the distance modulus they derive, /iNGC335i — 
30.01 ±0.19, is almost 0.4 mag closer than the Tanvir et al. 
(1995) distance for NGC 3368, MNGCases = 30.37 ± 0.12. 
This corresponds to a 2 Mpc difference in the line of sight 
distance to these galaxies, even though their projected sep- 
aration is 41 arcmin, or about 120 kpc at the inferred dis- 
tance. Graham et al. (1997) suggest one alternative ex- 
planation is that the /-band NGC 3368 Cepheid photom- 
etry of Tanvir et al. (1995) is systematically faint, which 
would lead to an underestimated extinction and an over- 
estimated distance. Metallicity corrections to the Cepheid 
distances do not change this conclusion, as both NGC 3368 
and NGC 3351 are similarly metal rich (Kennicutt et al. 
1998). If we adopt the mean My from the three other 
calibrators, SN 1998bu can in principle be used to test the 
Tanvir et al. (1995) distance. However, this procedure 
yields /iNGCsaes = 30.26 ± 0.27, and the uncertainty is too 
large to provide a definitive test. 

This discussion does point out an important attribute 
of our analysis. We have assumed that the internal errors 
in the Cepheid distances (i.e., excluding those factors we 
discuss in Section 5^) are accurately estimated. While we 
have tried to create a uniform set of distance moduli from 
the published values, a better procedure would be to re- 
analyze all the Cepheid data from the different groups in 
a completely consistent manner, beginning with the raw 
HST images. The Key Project team is doing this (Gibson 
et al. 1999); if needed, revised distance moduli for the SN 
calibrators can then be used in Table p^ to determine a new 
estimate for My and Ho- 

More information about the structure of the Leo I group 
can be determined by relying on the supernovae them- 
selves. For instance, the galaxy NGC 3389, host of the 
SN la 1967C, is sometimes considered a member of the 
Leo I group, though various group-finding algorithms dis- 
agree about its inclusion (Schneider 1989). We have used 
MLCS to analyze the photoelectric light curve of SN 1967C 
(de Vaucouleurs 1968), and conclude that NGC 3389 is at 
/x ~ 32.6, placing it significantly farther away than the 
other Leo I group galaxies. 

We can also compare the distance of NGC 3368 with 
that of NGC 3627, in the Leo triplet, located on the sky 



roughly 8 degrees from the Leo I group. NGC 3627 was 
the host of the SN la 1989B, which had a light curve sim- 
ilar to that of SN 1998bu. From our MLCS analysis of 
the SN 1989B we derive (FBrnax) = 12.01, Ay = 0.99 mag, 
and A = 0.19 mag, yielding m^ = 10.83±0.19. Combined 
with our result for SN 1998bu, m^ = 10.93 ± 0.18, we see 
that NGC 3627 is most hkely -- 0.1 mag closer than NGC 
3368, though there is a bit of uncertainty (arising from the 
correlated uncertainties in the MLCS analysis of each su- 
pernova). Since the supernovae only tell us the distances 
of their hosts, without further information we cannot de- 
termine the relationship between the Leo I Group and the 
Leo triplet as a whole. HST data has been taken to mea- 
sure a Cepheid distance to NGC 3627 directly; with our 
analysis of SN 1989B and taking M° = -19.34 ± 0.17 
derived from all four calibrators, we predict the distance 
modulus for NGC 3627, fi = 30.17± 0.25. In a very recent 
preprint, Saha et al. (1999) report their Cepheid distance 
to NGC 3627, with the result n = 30.22 ± 0.12, in ex- 
cellent agreement with our expectation. With this new 
Cepheid distance, we will be able to use SN 1989B as a 
calibrator in the future; including it in the analysis pre- 
sented here would decrease our estimate Ho by merely 0.3 
km s^-'^ Mpc""'^. 

This discussion strengthens our opinion that the sam- 
ple of calibrating SNe la should be restricted to those for 
which Cepheid distances have been directly determined, 
rather than relying on indirect association with other 
members of a group or cluster. Suntzeff et al. (1999) 
make this point with two SNe la in the Fornax cluster, SN 
1980N in NGC 1316 and SN 1992A in NGC 1380. Adopt- 
ing the Cepheid distance to NGC 1365 (Silbermann et al. 
1999; Madore et al. 1999) as the distance to the cluster 
and these two galaxies in particular leads to a correlated 
systematic effect on the SN la calibration. The analysis of 
Suntzeff et al. (1999) indicates that NGC 1365 is proba- 
bly foreground to both NGC 1316 and NGC 1380 by ~ 0.3 
mag. We have analyzed the SN la in both of these galaxies 
with MLCS and concur that NGC 1365 is closer than NGC 
1316 by - 0.3 mag. Our result for NGC 1380 is shghtly 
different, with the MLCS analysis suggesting NGC 1365 
is closer than NGC 1380 by only ^ 0.05 mag. Neverthe- 
less, this is exactly the sort of systematic error we wish 
to avoid, since the number of calibrators is small and sys- 
tematic errors in any of them can significantly affect the 
mean. 



6. CONCLUSION 

We have presented extensive photometric and spectro- 
scopic observations of SN 1998bu as well as an MLCS anal- 
ysis to determine the intrinsic luminosity (relative to the 
fiducial) of the supernova and the extinction along the line 
of sight. Using the Cepheid distance to NGC 3368 and 
three other SN la host galaxies, we have calibrated the 
absolute magnitude of our fiducial SN la, and applied this 
calibration to a set of 42 distant SNe la to derive a Hubble 
constant, Hq = 64^g km s""'^ Mpc~^, including systematic 
uncertainties such as in the distance to the LMC. The sta- 
tistical uncertainty in our estimate (~ 0.17 mag) arises 
from the fact that we have only 4 local calibrators; this 
uncertainty will be reduced by more SNe la occurring in 
galaxies with HST Cepheid distances, or more controUably 
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by measuring more Cepheid distances to the host galaxies 
of well-observed SNe la (SN 1998aq in NGC 3982 is an 
excellent target, for example). Reducing the systematic 
uncertainty in the Hubble constant (~ 0.15 mag) will be 
more difficult as it will entail a better understanding of 
the Cepheid distance scale and most importantly, a defini- 
tive distance modulus for the LMC. Nevertheless, we are 
optimistic that these reductions are possible, and eagerly 
await the day in the not-too-distant future when the Hub- 
ble constant is known to better than 10%. 
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